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This report oovers the seoond year of a two year study program, Initiated on 
July 9, 1979, conducted at Columbia University under NASA support, dealing with 
problems of multiple aooess, resource allocation and network design for next 
generation communication satellite systems. Two basic sets of requirements shaped 
our effort. First was the recognition that the next deeades will see increasing 
demand for the two limited resources at the disposal of the satellite system 
designer: suitable positions in geostationary orbit, and available radio frequency 
spectrum. Second was the fact that the user populations for these systems will 
be highly heterogeneous both in their statistical characteristics and their 
service requirements. Thus the general thrust of the work has been to attempt to 
identify the problems associated with the allocation of searoe communication 
resources in as efficient a manner as possible while satisfying the requirements 
of a heterogeneous and oontinually varying community of users. Particular 
attention has been directed toward systems involving large numbers of small 
"direct to user" earth stations. 

During the two year effort, we focussed on TDMA systems involving a high 
degree of frequency reuse by means of satellite switched multiple beams with 
varying degrees of onboard processing. The need for a system that oan provide 
integrated service to a disparate user group ranging from interactive computer 
terminals to television channels suggests that a highly sophisticated set of 
multiplexing, multiple aooess and switching procedures will be required if user 
needs are to be met economically while using the available bandwidth as 
efficiently as possible. Therefore, in the first year of the study we examined a 
variety of communication procedures, ohooslng those which showed promise of 
achieving these objectives, even though they might appear to be highly oomplex by 
today's operational standards. 
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The goal of the first year of the projeot was to develop and analyze the 
performance of multiple aooess techniques for satellite switohed multiple beam 
(SSf®) systems operating In a paoket switohed mode. The multiple aooess problem 
was viewed at two levels - the station aooess problem and the switoh assignment 
problem. Station aooess was oonsidered on a fixed (FA), demand (DA), and random 
access (RA) basis. For the switoh assignment problem, fixed (FS), demand (DS), 
and store-and-forward (SF) switching was considered. Six aooess protocols were 
selected for study, and delay- throughput analysis was performed for eaoh. One of 
the protoools, a demand-access/ demand-switched (DA/DS) "list scheduling" soheme, 
showed promise in a number of contexts, and was selected for further study in the 
second phase of our work. 

1.2. PROBLEMS ADDRESSED AND APPROACH USED 

As is typieal in any detailed analytical study, our first phase work raised 
at least as many questions as it answered. In the second year we addressed some 
of these outstanding issues. A major goal of this second phase was to examine 
performance questions in a more realistio traffic environment, using more 
realistic system configurations. In particular, we dealt with the following 
issues. 

1) Effect of Skewed Traffic, Overlapping Beams and Batched Arrivals in 

Packet-Switched SSMB Systems. 

2) Integration of Stream and Bursty Trafflo. 

3) Optimal Clrouit Scheduling in SSMB Systems . Performance Bounds and 

Computational Complexity. 
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The overall objective was to develop and evaluate algorithms for the 
effloient utilization of the satellite resouroes under reasonably realistio 
constraints. 

1.3 SUM4ARY OP RESULTS 

Simultaneous studies were undertaken to deal with the three problems 
mentioned above. The results of these three studies are reported in Chapters 2,3 
and 4 respectively. We summarize them in this seotion. 


1.3.1 Effects of Skewed Traffic, Beam Overlap, Batohed Arrivals 
Having observed from our first phase studies that the DA/DS "list scheduling" 
packet* switching protocol performed well in simple configurations, it was ohosen 
as a basis for further studies in more realistio environments. Three series of 
simulation studies were undertaken. (Due to the complexity of the system 
configurations analysis was out of the question.) A description of the protocol, 
the simulation approaoh and the techniques of validation are presented in Seotion 
2.1, with details and interpretation of results appearing in Seotion 2.2. The 
system parameters that were varied were: 

Number of zones - N 

z 

Traffic skew ratio - S 


Message length - I 


v. 


S, single packet 

M, mixed single arj multipacket 


Overlap probability 


f High 
i Low 




- 1.4 


In each ease the nuober of transponders was fixed at ■ 5 . 

The quantities of interest for perfomanoe evaluation were 
Average paoket queue length, Q 
Max queue length, Q njax 
Average paoket waiting time W 

m 

Average number of transponders used, H 
Max number of transponders used, N M . 

Each of these quantities was studied as a function of traffic intensity, c . 
The results of these runs are presented in a series of tables and figures in 
Section 2.2. 

The first study dealt with 1 small system (N z * a 5 ) , with no zone -verlap 
and single packet messages, the objective being to study the effect of traffic 
skew in isolation. 


For ease of reference, we list below the cases studied as well as the table 
and figure numbers where the results can be found. 


Case Mo. 


Skew ratio 
1 
2 


Table 

2.5 

2.5 


FlgSj. 


2.22-2.37 


3 


4 


2.6 


In the second study, we Increased to 20, keeping N^ * 5. In this way we 

were able to ooapare systems with N^ * (the aforementioned study) to those with 

"N z N t "* Basically, we were moving from a system which was "switch-limited” to 
one which was transponder- limited". Here, we considered the combined effects of 
skew and variable message lengths. 


Case No. 

Skew ratio 

Message length 

Table 



1 

5 

M 

2.7 

> 


2 

1 

M 

2.8 


2.41-2.47 

3 

5 

S 

2.9 


2.59-2.50 

4 

1 

S 

2.10 

J 



In the third study, beam overlap was considered in the context of a system 
with high traffic skew (S s 5 ), with N 2 s 20, N^ * 5. High or low overlap 
probability was created by skewing the traffic matrix more toward the overlapping 
or non-overlapping beams. (See Figs. 2.39—2.40) 


Case No. 

Overlap 

Message length 

Table 



Figs. 

1 

High 

M 

2.14 

\ 



2 

Low 

M 

2.13 


r 

2.48-2.58 

3 

High 

S 

2.12 



2 . 61 - 2.62 

4 

Low 

S 

2.11 

J 




As expected, the paoket queue lengths and average delay increased with 
trafflo skew, beam overlap and sassage length. However, because of the 
interaction of the sany parameters involved, the exact dependence of these 
quantities on the system configuration and traffic intensity was rather complex. 
One important result was that the combination of traffic skew and overlapping 
zones caused the system to saturate at lower traffic intensities than in systems 
with uniform traffic and no overlap. The most severe case studied (S = 5 with 

high overlap probability) resulted in p saturation 25 

1.3.2 Integration of Stream and Bursty Traffic 

In Chapter 3 we present an integrated circuit- switched and packet-switched 
multiple access scheme to handle combinations of stream and bursty traffic. The 
access problem is treated using a three-level hierarchical structure: level 1 is 

the concentration of the various information sources at their earth stations, 
level 2 is the multiplexing of the competing population of earth stations into an 
integrated TDMA frame, and level 3 is the allocation of satellite capacity and 
scheduling of stations to achieve some desired performance. 

The level 1 concentration problem has already been studied at length in the 
context of terrestrial networks. In this study we deal only with levels 2 and 3. 
An approximate analytic approach for predicting the behavior of the access scheme 
at level 2 is presented in Section 3.2, and it is used as a basis for performance 
analysis of the capacity assignment strategies discussed in Section 3.3. Four 
schemes for assigning circuit-switched (CS) slots are examined and compared: 


a) Fixed assignment: each station is allooated a fixed number of slots. 

b) Shared pool: some slots are allooated to the stations on a fixed basis, 
and others are assigned from a common pool on a demand basis. 

c) Demand assigned: all slots allocated on a demand basis from a common pool. 

d) Shared channels: a refinement of (b). 

Superimposed on the circuit- switched assignment is the packet-switched 
allocation. Ve stress here the movable boundary approach, in which the remainder 
of the frame (after CS capacity allocation) is reserved for packets. In addition, 
any CS 3lots not currently being used for active circuits are allocated to 
packets. 

In the integrated system, the parameters of interest are circuit blocking 
probability and packet time delay, which we examine and compare as a function of 
CS and PS traffic intensity. A numerical example is presented in Section 3*3*2 
which illustrates the comparative performance of the four CS allocation schemes. 
The best scheme from the point of view of blocking probability is (c). On the 
other hand, since (c) leads to the highest channel occupancy by circuit traffic, 
it gives the poorest packet delay performance . Thus, there is a trade-off 
between CS and PS performance. In all cases, the movable boundary scheme offers 
considerable enhancement of PS performance over fixed boundary schemes. 

In concluding this part of our study we also examined certain additional 
issues associated with on board processing and storage. These are discussed in 
Section 3.4, and include performance of downlinks with onboard buffering, and 
packet flow control based on CS slot oocupancy. 


1.3.3 Optimal Circuit Scheduling 

Focussing attention on a purely circuit-switched system we have attempted to 
draw some general conclusions concerning the feasibility of optimal scheduling in 
the face of adjacent beam interference, as well as other constraints caused by 
special system configurations. While this work had not been completed at the 
conclusion of the contract, we did obtain a number of interesting preliminary 
results. A complete account of this study, which was carried to conclusion after 
the termination of the NASA contract, appears in [GOP 82] . 

The preliminary results, presented in Chapter 4, deal with optimal switching 
algorithms for systems with single and multi-channel beams and with different 
uplink/downlink zone configurations. In the case of systems with adjacent beam 
interference, the scheduling problem is too complex to be solved optimally in an 
operational system. We have therefore devised some good heuristics for dealing 
with interference by a combination of proper switch scheduling and the use of two 
orthogonal polarizat ions . 


REFERENCE 

[GOP 82] Gopal, I.S., "Scheduling Algorithms for Multi-Beam Communication 
Satellites," Ph.D. Dissertation, Department of Electrical Engineeing, 
Columbia University, May, 1982. 
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2. DEMAND AS SIGNED/ DEMAND SWITCHED PACKET 

2.1 SIMULATION: DESCRIPTION AND VALIDATION 
2.1.1 Introduction 

This section addresses the effects of skewed vs. uniform traffic on 
packet delay of satellite-switched, time-division multiple-access (SS/TDMA) 
systems using the demand station access DA and demand switch assignment DS 
system as shown in Fig. 2. 1 . In the DA/DS system* all stations send request 
messages to the onboard scheduler for transmission slots for arriving data 
packets. Each request for slot allocation for the data includes information 
about the desired uplink-downlink connection. The scheduler makes the switch 
assignment on demand on a slot -by-slot basis, based on the list of requests 
currently awaiting scheduling, and triggers packet transmission from the 
stations whose requests have been scheduled. The operation of the DA/DS 
protocol is illustrated in the timing diagram of Fig. 2*2 • A data packet D^ 
arrives at a station which transmits a request on the order wire, indicating 
the desired destination. The scheduler forms a queue of all requests. The 
request waits on this queue for a time until it is scheduled. A 

signal T^ is transmitted by the scheduler on the downlink control channel 
to trigger transmission of the packet. On receiving the trigger the station 
transmits the packet. The average queuing delay is a function of the algorithm 
used by the scheduler for switch assignment. Simulation results are presented 
for the list scheduling algorithm briefly described in the next paragraph. 

The switch assignments must be made subject to the switching constraints which 
state that any one uplink zone may not be connected to more than one downlink 
zone at the same time. Similarly, any one downlink zone may not be connected 
to more than one uplink zone. Also, the total number of connections at any 
given time cannot exceed the number of transponders, N^. 

*See [STE 8(5} for further description of the DA/DS protocol, using the list 
assignment algorithm for processing requests. 


There can be many different algorithms for DS system, e.g. , maximum 
matching scheduling, cyclic permutations scheduling and list scheduling. 

An efficient algorithm is one which is simple to implement and which also 
reduces request waiting times. For our simulat ion we have considered the list 
scheduling algorithm, which is simple to implement though it ma not completely 
minimize the average waiting time in the request queue. We consider first a simple 
system in which the mmber of zones is equal to the number of transponders 
N t » 5. The contents of the queue of all pending requests at the controller 
is shown in Fig. 2.3 in matrix fo^m. Each cell in the matrix corresponds to 
a pair of source and destination zones; an entry in the cell indicates a 
request for a packet transmission between those zones. The total number of 
entries in the matrix indicates the length of the request queue, and the entry 
number indicates the position of the request in the queue. The switching 
constraints require that no more than one entry be scheduled in a given 
time slot from the same row or column of the matrix. In the list scheduling 
algorithm, the scheduler scans the lists of requests in the queue and schedules 
them if consistent with switching constraints. The scanning of requests is 
done on a first-come, first-served basis and terminates if either N^, requests 
have been scheduled or if all requests have been scanned. For the example of 
Fig. 2.3 a list scheduling of requests 1,2,8 is illustrated in Fig. 2,4 and 
Fig. 2 .5. 

We consider now a system with antenna arranged to form spot beams. 

Each beam comprises an uplink and a downlink at different carrier frequencies, 
with connections made between up and downlinks through a switch matrix in 
the satellite shown in Fig. 2.6. There are frequency translating trans- 
ponders, and each connection represents a path through one of them. Messages 


- 2.3 - 


arrive randomly at the stations and we assume that the message arrival processes 
at the stations are independent Poisson process. Each message consists of a 
single packet requiring one time slot for transmission. The spot beams can 
be overlapping or non-overlapping as shown in Fig. 2.7 and 2.8. In this section 
we consider systems with non-overlapping spot beams. 

Messages are assumed to arrive according to Poisson distribution with 
different mean interarrival time. Traffic can be uniform or non-uniform 
(skewed) for the zones of interest. The time slot length is assumed to 

be 400 units of time. 

2.1.2 Simulation Technique 

Discrete event simulation by digital computer offers many conveniences 
that make it an attractive tool for analysis of SS/TDMA systems. Simulation 
is essentially the technique of conducting sampling experiments on a model of 
the system. [kOB 78, DAN 81, FIS 78, SCH 74J . Thus simulated experiments 
may be regarded as similar to ordinary physical experiments and should be 
based on sound statistical techniques. Many simulation programs today are 
written in one of the existing simulation programming languages rather than 
in a general purpose programming language such as ALGOL, FORTRAN, PL/1 or 
APL. Some of the best known languages for discrete event system simulation 
are GPSS, SIMSCRIPT, SIMULA, CSL & SIMPL/1. Simulation language (1) provide 
a convenient representation of elements that couroonly appear in simulation 
models, (2) expedite configuration changes of the systems model, (3) provide 
an internal timing and control mechanism required to execute a simulation run, and (4) 
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facilitate collect ion of data and static ties on tha aggregate behavior of the I 

simulated system. In our simulation w have used the most frequently used 

language - GPSS. One of the special features of GPSS is that it describes 

directly the functional flow of the jobs (called transactions) through the 

system. The static structure of the system is described by entities such as 

facilities and storages. Facilities are entities that can be occupied by 

only one transaction at a time, while storages are entitles that may be 

occupied by more than one transaction simultaneously. 

Transactions have numerical characterlsitics called parameters. 

Facilities and storages have .lumerical attributes. These parameters and 
attributes are collectively termed standard numerical attributes. The dynamic 
mechanism of the system is modeled in terms of block commands: Each block 

command represents an operat ion that acts on individual transactions or 
entities. The operation is accomplished by means of field information or 
a collection of operands that accompany each block statement. 

Corresponding to each block command is a block-type subroutine. A GPSS 
simulation is the Interpretive execution of these subroutines. In order to 
execute simulated events in the correct time sequence, GPSS maintains a clock 
that records the current simulation clock time. All clock times in the 
simulator are integer values. The GPSS clock is updated to indicate the 
time of occurrence of the most imminent event. 

2.1.3 Simulation Model 

The basic flow chart for simulation is shown in Fig. 2.9. The model is 
composed of a major segment and two supporting segments shown in Figs. 2.10- 
2.12. These segments are described below. 


I 


Modal Segment 1 : Taking into consideration the most useful modeling approach, 
a transaction was used to simulate a packet. First we bring a packet into 
the model through use of the GENERATE block. The Interarrival tlsie is 
indirectly specified through halfword savevalue UTIL. This m a k es it possible 
to vary the traffic rate by redefining the value of UTIL between consecutive 
runs. Two byte parameters are associated with each transaction (packet), 
one assigned to a source zone and the other one to a destination zone. The 
function EXPON represents a continuous, 24 point GPSS function, used to sample 
from the exponential distribution with an expected value of 1. The function 
ADDR represents a discrete function for traffic distribution. The arriving 
packets form a queue. We use a QUEUE block in order to gather statistics. 
Since the time slot comes at regular discrete intervals of time and the 
packets can come at any time, we use a LOGIC GATE to let requests be scheduled 
only when a time slot arrives. The gate opens when a LOGIC SWITCH is set 
and the LOGIC SWITCH is set when a time slot arrives. Arrival of a time 
slot is shown in model segment 2. When a time slot arrives, we test for 
switch assignments of the requests in the queue. We use 11 Byte savevalues 
for this testing purpose. If a packet cannot be sent, it goes back into the 
queue with priority ahead of latter packet arrivals and waits for the next 
time slot. Since there are 5 transponders available, we use the INTER block 
to seize storage representing the transponder. 

Model Segment 2 : In this segment each transaction represents a time slot. 

As a time slot arrives, LOGIC SWITCH 1 is set, which causes the gate in model 
segment 1 to open. After a delay of 1 unit of time, the logic switch is 
reset. Since the assignment is done on a slot by slot basis, we initialize 
the Byte savevalues at the end of each time slot. 
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Model Segment 3 ; This Is a timer segment. Since the steady-state behavior 
is been considered here, statistics obtained during an initialisation 
period have been suppressed. The simulation is run for 2,000,000 units of 
time. 

Table 2.1 summarizes the definitions of GPSS entities used in the simulation. 


which is described further in Figs. 2.10-2.12 
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Table 2.1 DEFINITIONS FOR SIMULATION MODEL 


GPSS entity 

Interpretation 

Transactions 


Model Segment 1 

Packets 

Model Segment 2 

Time slot 

Model Segment 3 

A timer 

Storages 


1 

! 

Storage used to simulate switch connection 

Queues 


i 

i 

The queue used to gather statistics for 
the requests to be scheduled. 

i 

Functions 

EXPON 

Fucntion for sampling from exponential 
distribution with a mean value of 1, used 
to generate a Poisson arrival process. 

ADDR 

Function describing the distribution of 
traffic. 

Savevalues 

6 

Counter used to record the number of 
packets which have already been scheduled 
in a time slot. 

1-5, 7-11 

Used to store source zone and destination 
zone addresses. 

Variables 

UTIL 

Variable describing the mean interarrival 
time of packets. 

Logic Switches 

1 

A logic switch simulating the arrival of a 
time slot. 
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In order to achit-'e valid simulation reaules, it la naceaaary to aaaura 
that ataady-atata conditions ara approached and that the resulting statistics 
are statistically significant at a certain level of confidence. These con- 
siderations are treated below. 

2.1.4 Steady State 

A system has reached stable or steady state conditions when successive 
observations of the system performance are statistically indistinguishable. 

A system whose behavior does not satisfy steady-state conditions is usually 
described as being in a transient state. [ROB 78, DAN 81, FIS 78] For some 
kinds of systems no steady -state conditions are expected. 

Most simulation models are developed and run under the assumption that 
there are steady-state operating conditions. If this is the case, transients 
occur because the simulation run was begun with atypical values for its state 
variables. Therefore it is necessary first to select starting condition that 
keep the transient period short and then to minimise the effects of transients 
that do occur, so that the simulator performance is judged only on its steady- 
state behavior. 

For a single run of a simulator, a simple and basic strategy for setting 
starting conditions is to begin in an empty and idel status; that is, assume 
that the system is completely clean of activity and run the simulator until 
the transient effects are insignificant. It is easy to start the simulator 
under these conditions, but the transient period iu likely to be quite long. 

The effects of transients must be removed when performance measures 
are recorded in a simulation run. The most common method and the one used 
here is to introduce a nonrecoi ding interval to put the simulator into 
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steady-state conditions bsfors col lsc ting dsts. This Is dt«» by running ths 
slaulstor until steady-state conditions srs achieved . then clearing all 
statistical accunulatlons (but leaving ths stats of ths simulated systsn 
as it Is) and then continuing ths run. Ths conditions at ths end of ths 
transient period become an a priori estimate of ths steady-stats conditions 
and, in effect, are used to start a new run. This method is often easier 
to program than its alternative: to explicitly insert steedy-state starting 
conditions into the model. In GPSS, a RESET control operator facilitates 
the clearing and restarting process _GPS 71 ], This operator, which clears 
statistical mid other output accumulations without changing the model. Is 
used in our program. 

To determine whether the model has reached a steady-state condition is 
generally not a simple matter, because a steady state of a stochastic system 
or its model must be defined In terms of the probability distribution of 
the system state (rether than a particular state). That Is, the notion of 
the transient in a stochastic system is different, for Instance, from that 
of a deterministic system (suCh as an electrical circuit), where the transient 
behxvior is completely characterized by its impulse response function, it 
should also be recognized that the equilibrium Is a limiting condition that 
may be approached but never attained exactly. There is no single point in 
simulation time beyond which that system is in equilibrium, but we can choose 
some reasonable print beyond which we are willing to neglect the error that 
is made by considering the system to be In equilibrium. 

The choice of reasonable starting conditions can diminish the transient 
Interval but cannot eliminate it completely. No simple rules can be given 
to decide how long an Interval should be eliminated. It is advisable to 
use some pilot runs starting from the idle state to judge how long the initial 
bias remains. This can be done by plotting the measured statistic against 
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run lwgth u has beta dona In Figs. 2.13-2.15. Ic is highly desirable, however, 
that the pilot Investigation be dona by repeating runs. Needless to say, 
the initial values (often called "seeds") of the pseudo-random generator 
(used to simulate random arrival of packets and random source destination 
pairs) should be different In different runs; otherwise exactly the same 
system output would always be observed, Figs. 2.13-2.15 show how difficult it 
is to judge from a single run when the measured value has approached its 
steady state. To minimise simulation coat, not too many runs were made to 
determine steady stats. We have used only one random generator and since 
run length is proportional to computer cost, we have taken measurements for a 
certain limited interval of time only instead of making the run length too 
long. 


2.1.5 Confidence Intervals 

The use of confidence Intervals is an important technique for measuring 
the degree of uncertalnlty in making inferences about population character - _ -ice 
from sample data. The sample mean waiting tlma in queue W Is different from 
the actual population mean u and we estimate an interval around W such that 
we can predict with some level of confidence that u falls within that Interval. 
Out of the three main methods (KOB 78 for determining confidence intervals, we 
have used "the single run method" or "method of batch means". It makes 
consecutive runs to obtain the confidence intervals shown in Figs. 2.16-2. 18, 
using the final condition of one run as the steady-state starting condition 
for the next run. That is, one long continuous simulation run is divided into 
contiguous segments or batches, each of which has the seme length. The 
sample mean W^ of the j ctl segment is then treated as an individual obser- 
vation. The method of batch means is shown below. A simulation program was 
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written end the Initial 50,000 time units portion of output data starting form 
an empty state was discarded for the transient period. Pig. 2.19 shows the 
effect of simulation time on packet delay. A single run of T ■ 2 ,000,000 time 
units was divided into m*20 segments of equal length t a — • 100,000 units. 

Itl 

The observed process is the waiting time in queue W(t) as a function of 
simulated time t. More precisely, let the waiting time for a packet arriving 
at t ■ t^ during the j** 1 segment be denoted The sample mean of 

this quantity in the j 6 * 1 segment is then 

(j) 

- (j ) . 


where nj is the number of packets processed during the j segment , the 
overall sample mear. of the waiting time is 

m — (M 


( 2 . 1 . 2 ) 


and the sample variance is 


4 • > 2 


(2.1.3) 


The random variates W VJ 'are approximately normally distributed. Thus, the 
t-dlstribution of m-l"19 degrees of freedom [KOB 78, BAS 66] is our model for 
the variable 


U - u 


(2.1.4) 


where y is the (unknown) population mean waiting time E[w(t)] . 

The distribution of the variable t . is student's t -distribution with 

m-i 

m-i degrees of freedom. Fig. 2.20 shows the distribution curves for various 
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degrees of freedom. From these curves It can be seen thst with smell sample 
*1**, the distribution is considerably more spread out then the normal dis- 
tribution, but as 3ample size incresses the t-distrlbution approaches the 
normal distribution. By letting t Q ^ denote the upper j x 100 percentile 
of the t-distribution we have. 


-_ t < c y. - h> < t 1 - 1 - 

a/2;m-l — SL, — a/2;m-lj 


or 


P W - t , — — < p < W 

L o ' 2 '’*- 1 /T ” 


+ t 


a/2 ;o-l 


■ 1 - a 

(2.1.5) 

1 

— - 1-a 

(2.1.6) 

j 



The random variable W + c a /2* m -l ca ^ e< * a confidence interval and 

1- a is the corresponding confidence coefficient or the confidence level . 

The sample mean and variance based on the m(*20) data, uniform traffic 
with p ■ 0.2, are 


and 


W - 300-5 



m — {i ) — 2 

Z (W VJ '' - W) - 1127 
j-1 


(2.1.7) 


Thus, the confidence interval for this case is 


W + t „ . — — - 300-5 + 7.5 t 1Q . 

- a/2 ;m-l - a/2; 19 


( 2 . 1 . 8 ) 


Taking 0-95 as a confidence coefficient, for which tg Q25-19 “ 2-093, 
the 95 per cent (» 1 - a) confidence interval is calculated as 

w + t /<>.* i - JL - " 300-5 + 15.69 

- a/2 ;m-l ^ 


- 316-2, 284-8 

or 0.79, 0.712 slots 


(2,1.9) 


Similarly Cha 75 par cant confidence interval* for nonuniform traffic are 
also calculated. The results are given in Table 2.2, which shows that the 
mean waiting time obtained by simulation may be considered as reliable 
estimates. For the case with single packets only and uniform 

traffic. Table 2.3 and Fig. 2.21 shows the CPU time necessary to achieve 
reliable simulation results, as a function of traffic intensity. It will 
be noted that simulation of a large number of cases can become expensive. 
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Table 2.2 CONFIDENCE INTERVALS 


Traffic 

\Traffic intensity 
\ 0 

95% confidence 
interval (slots) 

Mean from simulation, 
W (slot 8) 

Uniform 

0.2 

0.71, 0.79 

0.76 

•1 

0.4 

1.12, 1.23 

1.19 

11 

0.6 

2.03, 2.25 

2.17 

Nonuniform 

(2:1) 0.2 

0.77, 0.87 

0.83 

•• 

(2:1) 0.4 

j 1.35, 1.52 

1.45 

it 

(2:1); 0.6 

3.24, 3.63 

3.49 

Nonuniform 

(4:1) 0.2 

0.82, 0.92 

t 

0.88 

tl 

(4:1): 0.4 

! 1.66, 1.94 

1.82 

t« 

(4:1) 0.6 

6.32, 7.63 

7.23 





■ . tv.> 

Table 2 .3 

UNIFORM TRAFFIC (Nj - N 2 - 5, single packets only) 


TRAFFIC INTENSITY! 
P 

CPU TIME IN 
MINS. 

i 

0.16 i 

0.39 

0.26 

0.56 

0.4 

0.85 

0.6 

| 1.50 

0.8 

3.04 
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2.2 EFFECTS OF SKEWED TRAFFIC, OVERLAPPING BEAMS, BATCHED PACKET ARRIVALS 
Tables 2.4. through 2.6 give the output statlstlos obtained for N^ • 
s 5 and Tables 2.7 through 2.14 give the output statlstlos obtained for N z 

= 20, N T » 5 from simulation. In general, various different kinds of 
traffics have been considered, 

i) uniform traffic 

ii) nonuniform traffic 

ill) nonoverlapping zone traffic 

iv) overlapping zone traffic 

v) single packet traffic 

vi) multiple packet traffic 

For nonuniform traffic, in the = N^. = 5 case, two types of traffic 
skew (2:1) and (4:1) have been considered, whereas for the - 20, N T * 5 case 
only one type of traffic skew (5:1) has been dealt with. With these 
different skews, the probability of source and destination zones are given in 
the corresponding tables. Nonoverlapping and overlapping zones for Nj * 20, 
Nt = 5 case are shown in Figs. 2.38 and 2.39, respectively. Two different 
kinds of overlapping zones have been considered. In Fig. 2.39(a) zones 1,2 
and 3 overlap as well as zones 6 and 7. This means that no two packets with 
source or destination address from the overlapping zones can be sent in the 
same time slot, e.g. , if zone 1 is sending or receiving a packet, zones 2 md 
3 cannot send or receive, respectively. On the other hand, if zone 2 is 
sending or receiving a packet, zone 3 can send or receive but not zone 1, 
since zone 2 overlaps with zone 1 and not zone 3- However, if zone 2 is 
sending, zone 1 may receive because the uplink and downlink frequencies are 
different. The probability distribution of the nonuniform traffic (5:1) for 
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N z s 20, N~, i 5 is shown in Figure 2.40. In Fig. 2.39(b) zones 13, 14 end 
15 overlap as wall as zones 18 and 19* The differenoe between Figs. 2.39(a) 
and 2.39(b) is that the traffio from overlapping zones in Fig. 2.39(a) is 
more probable than traffic from overlapping zones in Fig. 2.39(b). In the 
simulation results shown in various Tables and figures, traffio as in Figs. 
2.39(a) and 2.39(b) has been designated as overlapping zones (high 
probability) and overlapping zones (low probability) , respectively. Messages 
may consist of various numbers of packets. However, for simplicity, only two 
different kinds of message length statistics have been considered •messages 
consisting of only single packets and messages consisting of either single 
and or multiple packets each with equal probability of occurrence. In 
general, multiple packet messages might have many different lengths, but for 
our purposes, only messages with length eight times that of a single packet 
have been considered. In the simulation results we have represented messages 
with a combination of equiprobable single and multiple packets as "multiple 
packets." 

Tables 2.4 through 2.6 and 2.7 through 2.14 show the queue statistics 
and switoh matrix statistics in parts (a) and (b) respectively. The results 
of these various tables are plotted in Figs. 2.22 through 2.37 and in Figs. 
2.4l through 2.62. In any table for queue statistics, the most important 
parameters of interest, the average number of packets waiting in queue Q, 
mean packet waiting time in slots W and maximum number of packets waiting in 
queue have been shown as a function of traffio intensity p . Traffio 

skew s:1 Implies that the probability of having packets from a source 
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or to a destination is s times higher than that from other souroes or to 
other destinations. In the case of uniform traffio, skew is 1 : 1 . For the 
oase of N 2 > N t * 5 we have considered traffio skew 2:1 and 4:1; Traffio 
skew 2:1 means that Prob. [zone 1 to 3] *2 Prob. [zone 4 to 5] and skew 4:1 
means Prob. [zone 1 to 3 ] s 4 Prob. [zone 4 to 51 . For N z 3 20, N T * 5, 
traffic skew 5:1 has been considered, which means Prob. [zone 1 to 4] 35 
Prob. [zone 17 to 20], as shown in Fig. 2.40. 

The average number of paokets waiting in queue Q is given by the ratio 
of cumulative time integral to relative clock time [GPS 71] . The cumulative 
time integral for a queue is the sum of all of the clock units spent by all 
of the packets in the queue. The mean packet waiting time W is given as the 
cumulative time integral for the queue divided by total entries into the 
queue during the run. The maximum queue length 0 gives the maximum number 
of packets in queue at any time during the run. 

In any table for switch matrix statistics, the most important parameters 
of interest, maximum number of transponders used N max and average number of 
transponders used N, have been shown as a function, of traffic intensity 
P • N oax represents the maximum number of transponders used at any time 
during the run. The cumulative time integral for storage is the sum total of 
all of the clock units spent by all of the packets in storage. The average 
number of transponders N is given by the ratio of cumulative time integral 
to relative clock time. The traffic intensity p is defined as mean packet 
arrivals/slot/ transponder. 

A OA/DS system w:.th list scheduling has been considered. Now, 
according to the switching constraints, no more than one packet with same 
source (or same destination) zone can be transmitted in a given time slot. As 
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a result, in general, for nonuniform trafflo, the queue lengths and packet 
delays are expected to be larger than those for uniform trafflo. In oases 
of nonuniform trafflo, transmissions from the souroe or destination zones 
having higher probability mass functions are more likely to ooour, which 
implies that higher numbers of packets have repeating source (or 
destination) zones as ooopared to the oase of uniform trafflo. Due to our 
switching constraint, the larger the number of packets with the same source 
(or destination) zones, the more the queue builds up and and the average 
delay of the packets increases. Packets from high probability source or 
destination zones remain unprocessed with higher probability. Tables 2.4(a) 
to 2.6(a) and Tables 2.7(a) to 10(a) verify what has been discussed above. 
Comparing Tables 2. 4-2. 6, it is obvious that the more nonuniform the 
trafflo, the higher is the average queue content and the average delay. 

Tables 2.4 through 2.6 deal with s 5, and Tables 2.7 through 

2.10 deal with Ng s 20, N T * 5. Figs. 2.22 and 2.41 show a comparison 
of average number of waiting packets Q vs. traffic intensity for various 
traffic distributions. Fig. 2.22 deals with Ng * 5, Hj * 5 whereas Fig. 
2.41 deals with Ng * 20, ^ > 5, Figs 2.23, 2.42, 2.43, and 2.44 show 

comparisons of average waiting time V for packets vs. trafflo intensity p 
for various traffic distributions. Fig. 2.23 has ^ > 5 while the 
others have Ng s 20, N^= 5. denotes the maximum number of packets 

waiting in queue at any time during the simulation run, and since we have 
assumed Poisson arrival message statistics, we expect Q oax to increase with 
traffic intensity p. This is verified in Tables 2. 4-2. 6 and also in Tables 
2.7-2.10. Plots of maximum queue length Saax vs traffic Intensity p for 
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various traffio distributions ars shown in Pigs. 2.24-2.29 and Pig. 2.45. 

The first five figures deal with * 5 and the last one with H z * 20, 

N T a 5. The maxinun queue length oan have only discrete values since it 
represents a nunber of packets. There is not much difference in queue 
statistics for low values of traffio intensity , beoauae for low traffio 
intensity queue oooupanoy is low and the effect of the switching constraint 
is less. As a result it does not natter nuoh whether the traffio is uniform 
or nonuniform. But for higher values of traffio intensity there are many 
paoket arrivals and the switching constraints have a more significant 
effect, resulting in a stronger dependence of the results on different 
traffic distributions. 

For uniform and nonuniform traffic saturation ocours at different values 
of traffic intensity p. Por unifora traffio p & 1, for the 
nonuniform (2:1) case P s#t ^ 0.8 and for nonuniforn (4:1), p *at 
These are all for N z a N T * 5. 

Tables 2.4 -2.6 and Tables 2.7-2.10 also give comparisons of the 

average number of transponders used N and the maximum number of transponders 

used N for different traffic distributions. These results are plotted in 
sax 

Figs. 2.30-2.37 and Figs. 2.46-2.47. Like the naximun nunber of waiting 
packets in queue Q max , the maximum nunber of transponders used H miiy can have 
only discrete values. In general, we expect all the transponders to be busy 

when the traffic intensity is high for i 5 because for high traffic 

intensity, the chances of having waiting packets from all different source 
and destination zones are better, Which causes more transponders to be used. 
But for N z s 20, N T « 5, we expect the maxinun number of transponders used 
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H nax to be 5 for all valuta of trafflo intensity p baoausa in tha sans time 
slot, five paokets out of N z e 20 oan easily be found with nonrepeating 
source or destination address. But if Mj * 5, the ohanoes of having 5 
paokets with different sou roe or destination zone address are low. Thus 
stays constant at 5 in Fig. 2.47 and slowly rises to 5 in Figs. 2.34-2.37. 

Of course, the maximum number of transponders used oan never exoeed N^» the 
number of transponders available. From the plots we may oonolude that the 
distribution of traffio does not have mu oh influence on the maximum or 
average number of transponders used . The two results most strongly affected 
by traffio distribution are mean packet waiting time W and average queue 
length Q. 

Now considering traffic from overlapping and nonoverlapping zones, we 
expect traffio from overlapping zones to have larger mean waiting time and 
average queue length titan traffio from nonoverlapping zones, unless the 
overlapping zones have very low probability mass funetion. For traffio from 
the low probability regions, not many paokets from overlapping zones arrive 
and it does not matter whether the traffio is from overlapping or from 
nonoverlapping region. But if the traffic is from overlapping zones with 
high probability then we may expect several packets from the overlapping 
zones to contend for servloe. No two packets from overlapping zones oan be 
sent in the same time slot. In addition, every packet must also satisfy the 
condition that the overlapping zones are not sending (or receiving) any 
packet in the same time slot. The more constraints on the packet to be sent, 
the fewer the paokets that oan be sent in the same time slot. 

Tables 2.11-2.14 give the queue statistics for traffio from different 
kinds of overlapping zones and verify the conclusions above. These results 
are plotted in Figs. 2.48-2.52. From these simulation results it may be 
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concluded that for traffic from high probability overlapping zones, 
the seen packet waiting tine and the average queue o on tent is higher than 
that of low probability overlapping zones . However, by ooaparing Pigs. 2.42 
and 2.48 it may be concluded that the Man paoket waiting time for nonunifora 
traffio 3 : 1 ), nonoverlapping zones, and traffic from overlapping zones with 
low probability are very olose, especially for low traffio intensity. A 
similar conclusion about the average queue length may be drawn. 

Tables 2.11-2.14 also give comparisons of the average number of 
transponders used and the maximum number of transponders used vs. traffio 
Intensity for traffio from different kinds of overlapping zones. These 
results are plotted in Figs. 2.53*2.58. In these figures, the maximum 
number of transponders used does not depend on the type of overlapping zones. 
Sinoe M z « 20, N T * 5, there are always at least 5 packets satisfying all 
the constraints. The average number of transponders used is the same for 
different kinds of overlapping zones provided that all the traffio is served. 

However, saturation occurs at different values of p as oan be found 
from Tables 2.13 and 2.14. For multipaeket traffio from overlapping zones 
saturation ooours at P Mt v0.55 for high probability overlapping zones and 
at p 3at ££ .65 for low probability overlappi: * zones . 

Next the effect of message length on the queue of waiting paokets and on 
the transponder usage has been studied for * 20, a 5. Two types of 
message length statistics have been considered: only single paoket messages 

and a mix consisting of a combination of single and multiple paoket messages. 
Since acoording to the switching constraints, two packets with the same 
source or destination addresses oannot be soheduled in the same time slot, 
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we expect the queue of paoketa to build up for tha aultipla paokat 
oaaa and oauaa tha mean paokat waiting time W and tha average queue oontant 5 
to inoraaaa. Tablaa 2.7-2.14 giving tha quaua statlatioa for aingla and 
aultipla paokat traffic show this bahavior. Thasa results ara plotted in 
Figs. 2.42, 2.48, 2.59-2.62, whioh show near, paokat waiting tine W in slots 
for diffarant massage lengths and different traffic distribution, in Figs. 
2.41 and 2.51, whioh show tha effsot of aultipla paokats on tha avaraga 
number of paokats waiting for various traffic distributions, and Figs. 2.45 
and 2.52, which give the maximum number of packets waiting on the queue for 
different message mixes and traffic distributions. In general, the following 
conclusions may be drawn from the results; for multiple packets the average 
queue content and the mean paoket waiting time is auoh higher than that of 
single paoket trafflo whether the traffic is uniform, nonunlfora, overlapping 
or nonoverlapping. For the particular oase of aultipla packet and trafflo 
distributions considered, multiple paoket arrivals have auoh more influence 
on aean paoket waiting time and average queue oontant than nonuniforalty of 
trafflo. 

Tables 2.7-2.14 also give the maximum number and average numbs* of 
transponders used for messages of various paoket lengths and different 
trafflo distribution . These results are shown in Pigs. 2.46, 2.47, 2.53 - 
2.58. There is not auoh difference in average number of transponders used N 
for various paoket lengths and toe maximum number of transponders used N mv 
remains the same as before, since all packets are served below saturation. 
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TabU 2.4(a) QUEUE 
UNIFORM TRAFFIC 
F z - Hr - 5 

■ Prob. ^source son* ij . Prob. _ destination son* j] 
- j • j for 1 < i < 5 , 1 < j <5 


Traffic 

intensity 

P 

Average queue 
length, 

Q 

Mean packet 
waiting tine 
(slots), W 

Maximum queue 
length, Q m£X 

.07 

| 

.21 

.56 

4 

.2 

.77 

.76 

11 

.4 

; 2.37 

1.19 

16 

.6 

i 6.54 

2.17 

25 

.8 

i 

1 21.82 

f 

5.43 

61 


Tabl* 2.4(b) SWITCH MATRIX 


UNIFORM TRAFFIC 
N z - Nr - 5 


■ Piob .^source son* ij • Prob. ^destination son* jj 
•j'j for 1 < i < 5, 1<J <5 


I Traffic intensity, 

S 


Maximum no. of 

Average no. of trans- 

transponders used, 
N max 

ponders used, U 


.07 

.2 

.4 

.6 


3 

4 


.37 

1.01 

2 

3.01 
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Ta bl* 2 ■ S (.■») QUEUE 

NONUNIFORM TRAFFIC ; Traffic ikw 2:1 
N z - Nr - 5 
Traffic skew a:l 

Pj. • Prob . [iiource am# i] • Prob . [ daa C lnat ion tone j 


a- 

(3s + 2)' 

Us + 2) 2 
1 

(3s + 2) 2 


for 1 < i < 3, 1 < J < 3 

1 ^ i, J i it 4 ^ j , i i 5 
4 < i < 5, 4 < J < 5 



Table 2.5(b) SWITCH MATRIX 

NONUNIFORM TRAFFIC; Traffic skew 2:1 

N - N - 5 
Z T 


Traffic intensity. 

Maximum no. of trana- 

Average no. of trans- 

P 

ponders used, N m>x 

' ponders used, N 


.6 

.63 


3.15 

3.34 

V6S 
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Table 2.6(a) QUEUE 


NONUNIFORM TRAFFIC; Traffic skew 4:1 


n t" 5 


Traffic skew s:l 

■ Prob. ^source zone l"* 


Prob, ^destination zone 


(3s + 2) 2 

3 

(3s + 2) 2 
1 

(3s + 2) 2 


for 1 <_ i ^ 3, 1 j je 3 

111, j < 3, 4 < j, 

4 < i < 5, 4 < j £ 5 


J 


1 < 5 


I 1 

Traffic intensity, ! 

P ! 

Average queue 
length, Q 

t 

Mean packet j 

waiting time 
(slots) , W 

| 

Maximum queue 
length, ^ 

.07 

.23 

.61 

5 

.2 

.89 

.88 

11 

.4 

3.66 

1.83 

19 

.6 

21.8 

7.23 

60 

.612 

! 

26.58 

8.62 

67 , 


Table 2.6(b) 
NONUNIFORM TRAFFIC; 


N„ 


N„ 


SWITCH MATRIX 
Traffic skew 4:1 


Traffic intensity, 

p 

Maximum no. of trans- 
ponders used, N_ 

max 

Average no. ofjtrans- 
ponders used, N 

.07 

2 

.37 

.2 

3 

1.01 

.4 

4 

2 

.6 

5 

3.01 

.612 

5 

3.08 
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Table 2, 7 fa) queue 

NONOVERLAPPING ZONES 
NONUNIPORN TRAFFIC; Traffic skew 5:1 
N z * 20 , N T « 5 

Pjj s P rob. [source zone i]* Probl destine tion zone Jl 
Prob.[ source zone i] * Prob.[ destination zone J] 
s 0.083 for 1 < i, j < 4 
» 0.066 for 5 $ i, j £ 8 
s 0.049 for 9 < i, J < 12 
3 0.033 for 13 < i, j < 16 
s 0.016 for 17 £ i, J £ 20 


MULTIPLE PACKET 


Traffic intensity, 

p 

Average queue 
length, S’ 

Mean packet waiting 
time (slots), W 

Maximum 
queue length, 

®max 

0.2 

3.88 

4.02 

32 

0.4 

9.22 

4.57 

43 

0.6 

16.10 

5.40 

67 

0.8 

30.50 

7.59 

108 


Table 2.7(b) SWITCH MATRIX 
NONOVERLAPPING ZONES 
NONUNIFORM TRAFFIC; Traffic skew 5:1 


N z s 20, N T s 5 
MULTIPLE PACKET 


Traffic intensity, 

p 

Maximum no. of 
transponders used, 

"max 

Average no. of transponders 
used, N. 

i 

0.2 

5 

0.97 

0.4 

5 

2.01 

0.6 

5 

2.97 

0.8 

5 

4,01 

















- 2.28 - 


Table 2.8(a) QUEUE 
NONOVERLAPPING ZONES 
•jNIVORM TRAFFIC 
N z » 20, N T = 5 

Pjj = Prob-[ source zone i]. Probl destination zone j) 

ProfcJ source zone i] = Probi destination zone J] 

= __L. for 1 I i, j 1 20 
20 

MULTIPLE PACKET 


Traffic intensity, 

p 

Averag queue 
length, Q 

Mean packet waiting 
time (slots), W 

Maximum 
queue length, 

®max 

0.2 

.n 

3.92 

32 

0.4 

8.79 

4.36 

46 

0.6 

14.87 

4.99 

59 

0.8 

27.30 

6.80 

107 


Table 2. g(b) SWITCH MATRIX 
NONOVERLAPPING ZONES 
UNIFORM TRAFFIC 



MULTIPLE PACKET 


Traffic intensity, 

p 

Maximum no. of transponders 
UMd ’ "max 

Average no. oft rans- 
ponders used N 

0.2 

5 

0.97 

0.4 

5 

2.02 

0.6 

5 

2.97 

0.8 

5 

4.00 
















Table 2.9(a) QUEUE 


NONOVERLAPPING ZONES 
NONUNIFORM TRAFFIC; Traffio skew 5:1 


H z* 

jj * Prob. [source 
Prob.[ source zone 

20, N t *5 

zone il . Prob. (destination zone Jl 
i] s Prob. (destination zone j] 

= 0.083 

for 

1 < i,J i 4 

= 0.066 

for 

5 i i.J <_ 8 

= 0.049 

for 

9 5. i,J 1 12 

= 0.033 

for 

13 i i, j 1 16 

= 0.016 

for 

17 i i,j < 20 


SINGLE PACKET 


Traffic intensity, 

p 

Averse queue 
length, $ 

Mean packet waiting 
tine (slots), W 

Maximum 
queue length, 

®max 

0.2 

0.59 

0.58 

7 

0.4 

1.32 

0.66 

12 

0.6 

2.37 

0.78 

16 

0.8 

4.38 

1.09 

20 


Table 2.9(b) SWITCH MATRIX 
NONOVERLAPPING ZONES 
NONUNIFORM TRAFFIC; Traffic skew (5:1) 
N z * 20, Nj *5 
SINGLE PACKET 


Traffic intensity, 

p 

Maximum no. of transpon- 
ders used , N 

max 

Average no. of transpon- 
ders used, N 

0.2 

5 

1.00 

0.4 

5 

2.00 

0.6 

5 

3.01 

0.8 

5 

4.02 


















Table 2.10(a) QUEUE 
NONOVERLAPP IHQ ZONES 


UNIFORM TRAFFIC 
N z s 20, N t * 5 

P^j a Prob. [source zone il . Prob(destination zone jl 

ProtaJ source zone il » ProbI destination zone j] 

a_l_ for 1 < i,J <20 

20 “ 

SINGLE PACKET 


Traffic intensity, 

p 

Average queue 
length, Q 

Mean packet waiting 
time (slots), W 

Maximum queue 

0.2 

0.56 

0.55 

7 

0.4 

1.23 

0.61 

12 

0.6 

2.16 

0.72 

17 

0.8 

4.05 

1.01 

20 


Table 2.10(b) SWITCH MATRIX 
NONOVERLAPPING ZONES 
UNIFORM TRAFFIC 

N z a 20, N t * 5 

SINGLE PACKET 


Traffic intenstiy, 

0 

Maximum no. of 
transponders used, 

"max 

Average no. of transponders 
used, N 

0.2 

5 

1.00 

0.4 

5 

2.00 

0.6 

5 

3.01 

0.8 

5 

4.02 
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Table 2. 11 (a) QUEUE 

OVERLAPPING ZONES (overlapping in lot# probability region) 
NONUNIFORM TRAFFIC; Traffic skew 5:1 
N z * 20 , N T = 5 

~ * >ro * ) ‘ (source zone 1]. Prob. (destination zone Jl 

Prob. [source zone 1] s Prob. [destination zone j] 


= 0.083 

for 

1 

i i.J 

l 1 * 

= 0.066 

for 

5 

1 

1 8 

S 0.049 

for 

9 

i 1;J 

1 12 

= 0.033 

for 

13 

1 ij 

1 16 

= 0.016 

for 

17 

1 i,j 

1 20 


SINGLE PACKET 


Traffic intensity, 

p 

Average queue 
length, Q 

Mean packet 
waiting time (slots), 

W 

Maximum queue 
length, 0^ 

0.2 

0.59 

0.58 

7 

0.4 

1.33 

0.66 

12 

0.6 

2.43 

0.8 

16 

0.8 

4.57 

1.14 

21 


Table 2.11(b) SNITCH MATRIX 
OVERLAPPING ZONES (overlapping in low probability region) 
NONUNIFORM TRAFFIC; Traffic skew 5:1 
N z * 20, N T * 5 

SINGLE PACKET 


Traffic intensity, 

p 

Maximum no. of 
transponders used, 

N max 

Average no. of 
transponders used, 

N 

0.2 

5 

1.00 

0.4 

5 

2.00 

0.6 

5 

3.01 

0.8 

5 

4.02 
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Table 2.12 (a) QUEUE 

OVERLAPPING ZONES (overlapping in high probability region) 
NONUNIFORM TRAFFIC; Traffio skew 5:1 
N z a 20, N T * 5 

Pjj a Prob. [source zone i). Probi destination zone Jl 

Protxl source zone i] a Prob. [destination zone j] 


0.083 

for 

1 

< iJit 

0.066 

for 

5 

i ij £ 8 

0.049 

for 

9 

i 1J £ 12 

0.033 

for 

13 

I i,J £ 16 

0.016 

for 

17 

£ i, J 1 20 


SINGLE PACKET 


Traffic intensity, 

p 

Average gueue 
length, Q 

Meand packet 
waiting time (slots), 
W 

Maximum queue 
queue length, 

®max 

0.2 


0.65 

9 

0.4 

1.59 

0 (9 

14 

0.6 

3.26 

1.08 

19 

0.8 

7.32 

1.82 

'UL 


Table 2.12(b) SWITCH MATRIX 


OVERLAPPING ZONES (overlapping in high probability region) 
NONUNIFORM TRAFFIC; Traffic skew 5:1 


N z a 20, N T a 5 

SINGLE PACKET 


Traffic intensity, 

P 

Maximum no. of 
transponders used, 

N max 

Average no. of 
transponder used, 

N 

0.2 

5 

1.00 

0.4 

5 

2.00 

0.6 

5 

3.01 

0.8 

5 

4.02 
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Table 2.13 (a) QUEUE 

OVERLAPPING ZONES ( overlapping in low probability ragion) 
NONUNIFORM TRAFFIC; Traffic skew 5:1 
N z a 20, N T * 5 

Pjj s Prob. [source zone il . Prob. [destination zone jl 
ProbI source zone i] * Prob. [destination zone j] 


a 0.083 

for 

1 

I M 

I * 

a 0.066 

for 

5 

1 i ,3 

1 8 

a 0.049 

for 

9 

1 i,J 

<12 

= 0.033 

for 

13 

1 i» J 

1 16 

a 0.016 

for 

17 

I i,J 

120 


MULTIPLE PACKET 


Traffic intensity, 

p 

Average Queue 
length, Q 

Mean packet 
waiting time 
(slots), 9 

Maximum 
queue length, 

Saax 

0.2 

3.9 

4.04 

32 

0.4 

5.4 

4.67 

45 

0.6 

16.66 

5.60 

70 

0.65 

19.69 

6.02 

78 


Table 2 . 13 (b) SWITCH MATRIX 


OVERLAPPING ZONES (overlapping in low probabililty region) 
NONUNIFORM TRAFFIC; Traffio skew 5:1 


N z = 20, N T = 5 

MULTIPLE PACKET 


Traffic intensity, 

p 

Maximum no. of 
transponders used, 

N max 

Average no. of trans- 
ponders used , N 

0.2 

5 

0.96 

0.4 ; 

5 

2.01 

0.6 

5 

2.97 

0.65 

5 

3.27 









2.34 - 


Table 2.14 (a) QUEUE 

OVERLAPPING ZONES (overlapping in high probability region) 
NONUNIPORM TRAFFIC; Traffio skew 5:1 
N z » 20, N t * 5 

P iJ * Ppob *l* ouree zone U • Prob. [destination zone j] 

Probi source zone i] * Prob. (destination zone J] 
s 0.083 for 1 <_ i, J <_ 4 

* 0.066 for 5 li,J 1 8 

a 0.049 for 9 ii,J 1 12 
a 0.033 for 13 < i J 1 16 

a 0.016 for 17 i i,J 120 


MULTIPLE PACKET 


Traffic intensity, 
P 

Average queue 
length, 5 

Mean packet 
waiting time (slots), 

ff 

Maximum 
queue length, 

®max 

0.2 


4.3 

34 

0.4 

' ' EM* J v 

5.4 

59 

0.45 


5.65 

60 

0.51 

16 

6.21 

64 

0.55 

18.24 

6.58 

70 


Table 2.14(b) SWITCH MATRIX 


OVERLAPPING ZONES (overlapping in high probability region) 
NONUNIFORM TRAFFIC; Traffio skew 5:1 

N z a 20, N T a 5 

MULIPLE PACKET 



Maximum no. of 
transponders used, 

N max 

Average no. of 
transponders used, 
N 

0.2 

5 

oT?? 

0.4 

5 

2.01 

0.45 

5 

2.23 

0.51 

5 

2.57 

0.55 

5 

2.77 
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FIG. 2.2 PACKET SCHEDULING IN A DA/DS SYSTEM 
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Fig. 2.10 Model Segment 1 
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Fig, 2.13 Effect of simulation time on packet delay for uniform traffic* 
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Fiq. 2.14 Effect of simulation time on packet delay for nonunlform traffic (2:1). 
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Fig. 2.18 Confidence interval on waiting tijae for nonuniforn traffic (4:1) 
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Fig. 2.19 Effect of stipulation time on packet delay with initial 50,000 
units of simulation time as the nonrecording period. 






Fia. 2.21 Cost vs traffic intensity curve. 
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Fig. 2.24 Maximum number of waiting packets ((^ ) for uniform traffic. 
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Fig 2.25 Maximum number of waiting packets for uniform traffic. 
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Fie. 2.27 Maximum number of waiting packets (Q^) for ponunlform traffic (2:1) 
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Fig. 2.28 Maximum number of waiting packets ((^ ) for nonuniform traffic (4:1). 
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Fig. 2.30 Average no. of transponders used for uniform traffic. 
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Fig. 1.33 Average number of transponders used for nonuniform traffic (4:1), 
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Fig. 2.34 Maximum number of transponders used for uniform traffic. 
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Fig. 2.36 Maximum number of transponders used for nonuniform traffic (4:1). 
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Fig. 2.44 Comparison of waiting time for different traffic distributions. 
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Fig. 2.61 Comparison of waiting time for different message lengths. 
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3. INTEGRATION OF STREAM AND BURSTY TRAFFIC 
3J. INTRODUCTION 

3.U Conventional Segregated-Switched Networks: 

It has been a tradition in both public and private sectors to have 
communication networks handling separately voice and data traffic. 

Initially, the transmission of information was successfully implemented in 
two modes: Analog for voice, digital for data. Correspondingly, two 
switching techniques were developed and proven successful up to now: 
o Circuit Switching for Voice 
o Packet Switching for Data 

The circuit switching (CS) mode consists of establishing a source-to- 
destination physical path (circuit) and maintaining it for the entire duration 
of the transaction (call). The cross network delay is then small and 
approximately constant. The network manages the access of CS traffic on a 
loss basis; i.e. if the outgoing trunk is fully busy, then arriving calls 
(customers) are blocked. The degree of service in such networks is typically 
measured by the probability of blocking (or loss) of a given call . The 
packet switching (PS) mode consists of grouping data into packets, each 
packet being routed independently through the network on a store-and- 
forward basis. The network communication resources (Buffers, links, etc.) 
are dynamically shared, hence, better utilized. The degree of service 
expected from a PS network is commonly evaluated by two interrelated 
measures: The end-to-end delay suffered by a transmitted packet and the 
delivery rate of packets to their destinations, referred to as throughput. 

In this mode of switching, time delay and throughput are subject to severe 
degradation when the communication resources of the network are 
overloaded by an uncontrolled offered traffic load. To prevent or alleviate 
such per formanoe degradation and maintain the network operating point in a 


desired range of efficiency (minimum delay, maximum throughput), three 
issues have been actively studied in recent years, namely: 
o Flow control 
o Routing 

o Buffer management 

These issues, although important, are not addressed in this work, with the 
exception of section (3.4.1.2), where a simple data flow control mechanism is 
presented. 


3.1.2 Evolving Integrated-Swit chad Networks 
As the variety of users is rapidly expanding and the volume of their 
traffic is expected to increase Sharply, future communication networks are 
planned to have the capabilities to handle heterogeneous leers with different 
characteristics and requirements. Considering t/pical classes of traffic, and 
based on the transaction size, traffic rates and delivery time requirements, 
we may classify this traffic mix into two categories: 

o Stream Traffic: (e.g., voice, video, file transfer, ... ) 
o Bursty traffic: (e.g., Interactive computing, information 
retrieval »•* ) 

Stream traffic messages are mainly long and uninterrupted. Their 
users have a wide spectrum of characteristics (e.g., the bit rate ranges from 
tens of Kbps for voice to Mbps for video) and requirements (e.g., the max 
time delay variation ranges from 200 ms for voice to possibly tens of minutes 
for file transfer). This traffic may be managed on a blocking basis (e.g., 
voice) or on a delay basis (e.g., data). 

Bursty traffic messages are short and discrete with a high peak-to- 
average traffic ratio. The transmitted information requires extensive error 
control with tolerable total network delay. This traffic is commonly 
managed on a delay basis. 

In a communication network, the resource of primary cost is the 
transmission bandwidth available (or trunk capacity). The switching mode 
(CS or PS) for these 2 classes of traffic is thus based on the efficiency 
achieved in bandwidth utilization. Thus, CS is clearly more cost-effective 
for stream-type traffic because of the long and contiguous nature of 
information messages. CS leads to bandwidth waste during the inherently 
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long interruptions of bursty traffic flow. Hence, PS is a more likely 
candidate for bursty-type traffic. Figure 3.1 shows how CS and PS are 
compared in terms of time delay with respect to message length [OCC77] . 

As was mentioned earlier, CS and PS networks evolved traditionally in 
parallel to support separately voice and data communication needs, 
respectively. As more and more communication services emerge (FAX, video 
teleconferencing, word processing etc.), it is expected that the future trend 
will favor the integration of both CS and PS in a common rather than 
separate system. This integrated switching approach promises the following 
benefits (MAG 81] 

o Dynamic sharing of the network facilities leads to much higher 
utilization of communication resources, yielding significant dollar 
savings. 

o New technology trends and new applications are more easily 
adapted to integrated networks. 

o Equipment duplication is minimized because there is more 
flexibility in interconnections between the various types of user 
terminals. 

However, in contrast to these advantages, drawbacks may be imposed 
at the implementation level of an integrated system. This is due primarily to 
the large diff erenees between applications requirements (e.g., bandwidth 
ranging from tens of Hz for TTY to several Mhz for video, digital 
transmission for some applications, analog transmission for others...) 
Implementation costs, in both software and hardware, may be considerably 
higher and will have to be studied, however. 
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3.1.3 Overview of the Existing Integrated-switching Approaches 

Apparently, the earliest proposal for intagrated switching was initiated 
by IBM in 1974 [JAN 80) . It consisted of multiplexing CS and PS 
traffic in a slotted frame structure representing the common channel. 

Figure 3.2 illustrates a time division integrated multiplexor. The time 
axis is quantized by frames subdivided into time slots. A boundary separates 
the frame into a CS portion which handles CS traffic (e.g. , voice) and a 
PS portion which handles PS traffic. This is referred to as fixed frame 
fixed boundary (FFFB). Due to the unpredictable variations in CS traffic, 
part of the CS region of the frame may be unused, thereby wasting band- 
width. To avoid this wast *. and improve channel utilization, PS traffic 
may be allowed to cross the nominal boundary and occupy idle CS slots. 

This is referred to as fixed frame movable boundary (FFMB). A variable 
frame which adapts its length to traffic variations has also been considered 
in the literature. 

The first analytical investigation of these multiplexing structures 
(FFFB, FFMB) was carried out by Kum merle [JAN 80] who derived 
approximate 

expressions for circuit blocking probability (using classical Erlang-B-formula) 
and average packet delay. His results, supported by simulation, 
showed higher performance of the MB, as may be expected. In 1976, Fisher & 
Harris [FIS76] attempted an exact analytical approach to the MB system, 
using a multiserver model technique, resulting in an M/D/N queue 
formulation of the PS process (data). This analysis was validated by 
assumptions which were reported to be erroneous in [JAN80,WEI80] and the 
model is claimed to give very underestimated values for packet delay when 
the effective capacity available for PS traffic is smaller than the offered 
load. 
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In 1977, an approximate modal was developed by Occhiogroeso et al 
[OCC 771 , in which an arbitrary slot within a frame is intermittently 
available for packet service. Analysis of this model is tractable, and is 
presented in Section 3.2 for subsequent use. At about the same time, Chang 
approached the analysis of the FFMB system using an exact for- 
mulation of the joint voice/data stochastic process as a two-dimensional 
Markov chain. The analysis was successfully carried out and results were 
succinctly obtained, but the computational complexity poses a severe limi- 
tation in treating realistic cases. In order for this report to be self- 
contained, this approach is also presented in Section 3.2 and its complexity 
is pointed out. 

In the work of Maglaris Sc Schwartz [MAG 79] , the variable frame 
scheme is first analyzed and its superior performance is shown as compared 
to the fixed frame. However, its implementation needs considerable 
hardware and software complexity to meet synchronization and time- 
transparency requirements for CS services. The fixed frame multiplexing 
structure was modeled as a Markovian decision problem under two optimality 
criteria involving blocking probability for CS traffic and time delay for PS 
traffic. Again, according to this model, the MB scheme offers near optimal 
performance. 

Recently, Janakiraman [JAN 80] modelled the fixed and variable 
frame schemes both with fixed and moving boundaries, with a single server 
(synchronous switch) and gave a three-dimensional Markov chain formulation. 
The three dimensions are associated respectively with the CS process, PS 
process, and position of the server (switch) within a frame. Numerical 
results of this contribution agree quite with previous work. Other contri- 
butions in this area, mostly numerical (simulation), are also known in 
the literature. Table 3.1 shows a chronological list of the main contributions 
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with a summary of thair characteristics [JAN 80) . 

3.1.4 Problem Definition and Outline of the Chapter 
The 30/20 QHz market assessment studies prediet a large number 
of services included in three main categories! Voice, Data, Video. This 
heterogeneous population of users ranging from interactive computer termi- 
nals to television channels suggests a high integration of multiplexing, 
multiple access and switching procedures in order to meet economically the 
user needs with efficient usage of the available bandwidth. In this chapter 
we address the isues of integrated swiehing and bandwidth allocation in the 
environment of the 30/20 GHz multiple-beam satellite communication 
system. The problem of user access to this satellite system can be seen in a 
three-level control hierarchy. 

Level 1. Connection and Concentration of the various types of 
information sources (classified in our terminology into stream 8c 
bursty terminals) to their earth-station. 

Level 2. The M stations, competing for the same uplink channel, have 
to be multiplexed in an integrated structure. The TDMA protocol at 
this level is considered throughout this chapter. 

Level 3. The onboard switching facilities are shared by the 
zones (covered by a multibeam antenna) to access the Nz downlinks 
via Ny, transponders, Nj < Ng. 

A strategy is needed to sequence and schedule the 

connections of uplinks to the appropriate downlinks in order to achieve 


some desired criterion. 

This hierarchy is illustrated in Fig. 3.3. 
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The access problem of terminals to their earth-station (Level 1) is not 
addressed here. We consider first the integrated multiplexing of stations to 
access their common uplink channel (Level 2), second, we examine the 
problem of scheduling the uplink-downlink transactions, with attention paid 
to the (constraint of time requirement inherent to some CS traffic. 

In Section 3.2 the analysis problem associated with CS and PS 
stochastic processes arising in integrated systems is posed; then the queuing 
problem is analyzed by adapting two models that have appeared in the 
literature. These models have been selected to point out, on one hand, the 
computational intractability of an exact analysis, and on the other hand, the 
simplicity of an approximate method. 

In Section 3.3, we approach the problem of allocating satellite 
capacity to integrated users (stations) communicating via the satellite in a 
multizone, multibeam system. Focussing on the station access problem 
(Level 2), different capacity assignment strategies are presented and 
analyzed for CS traffic. A numerical example is treated and performance 
evaluation is carried out for both CS and PS traffic. The closed form 
results obtained from the approximate models of Section 3.2 are used for 
packet performance evaluation. Finally, the different capacity assignment 
schemes are compared in terms of both circuit and packet performance. 
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OF POOR QUALITY 

3.2 TRAFFIC BEHAVIOR AT THE INTEGRATED MULTIPLEXOR LEVEL 

We adopt the common assumption of Poissonian statistics for 
both circuit and packet arrival processes. In addition, with 
the assumption that both circuit switched call durations (circuit 
holding times) ai.J packet lengths (packet transmission times) 
are exponentially distributed random variables, we achieve also 
Poissonian statistics for circuit disconnection and packet 
departure processes. Thus, the joint circuit/packet switched 
(rp/ps) traffic process arising in an integrated multiplexor 
evolves at steady state as a two-dimensional Markovian birth- 
death process. The parameters of this model are described in the 
following section. 
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3.2.1 Analytic Model with Exponential Message Lengths: 
Two*dimensional Markov Chain 

The integrated system, described earlier, is essentially a 
queuing system with two types of arrival processes: stream arrival 

(CS) and bursty arrival (PS) . An FFMB frame management is adeoted 
in this analysis because of its higher performance (proven in 
relevant studies: FIS 76, MAC 79 etc) and its potential feasibility 
for use in an integrated-network design procedure. The integrated 
system in this case is succinctly characterised by the following 
characteristics and parameter assumptions: 

. CS stream traffic 

Arrival process is assumed Poissonian with rate 
Holding time is exponentially distributed with mean 
Managed on a blocking basis. 

. PS bursty traffic 

Arrival process Poissonian with rate *2* 

Packet length exponentially distributed with mean - . 

U2 

Managed on a buffering basis. 

Buffer size finite / infinite 

. Channel characteristics 

Frame duration: F slots 
CS portion : N g slots 

PS portion : slots 


<= |H* 
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. Performance measures 

CS traffic: Blocking probability Pg . 

PS traffic: Average pkt delay and/or probability 

of overflow. 

The time slots (i.e. circuits) are assigned to CS users 
(customers) at the start of each frame (SOD. It is assumed that 
CS customers arriving during a frame are queued until the start 
of the next frame and, if all the N,. slots are busy then these 

O 

customers are lost. In queuing terminology, the CS process behaves 
like an M/M 'N,. /\' s process. Packets are sorted in each of the Ng 
slots plus anv idle slot in the CS region. An arriving CS customer 
which cannot find a fre*> slot preempts a PS customer. Preempted PS 
customers return to the data butter and await service. The PS process 
behaves like an M/M 'Ng*S/M ? recess, where S is a random variable rep- 
resent it. g the number of idle circuit slots (0-S-N^.l and B is the buffer 
sire. This svstem, illustrated in Pig. A. I, lends itself to a two- 
dimensional Markov chain formulation, as in Fig. A . 2 and Appendix A. 

The variables and Qg represent the number of CS and PS customers in 
the svstem, respectively. The state diagram of this joint process and 
the corresponding state equations for the case B finite appear in 
Appendix A. Tn principle, the state equations could he solved to 
determine the steady-state joint probabilities for various numbers of 
occupied stream traffic (CS) slots anJ bursty (PS) buffer occupants. 
However, it is not practical to do so. (See Appendix A). 

Case B * » (Infinite buffer si:e) 

If the storage space available (B buffer units) for data is 
assumed to he verv large compared to the date queues that mav build 
up, we can consider the buffer sire (B) to he infinite. In this case, 
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the system of difference equations appearing in Appendix A, is solved 
by the use of the standard technique of moment generating functions 
(mgf). Given that taere are Qg * i CS customers in the system, we 
define the coiditional m.g.f. as: 

Pi(z) - fj = 0 for i=0,l,...N s (3.2.1) 

Applying this m.g.f. to the equations (A.l) to (A. 6) in Appendix A, we 
get the linear system 

p 0 m ' b n (z) 

V--) b i (z) 

Aii = (3.2.2) 

i 

! P NS (z) | b NS (z) ' 

I . s 1 

■ ^ \ 

Where Ai j , the (i,j)the element of the NgXNg system matrix, is a 

function solely of the parameters (A^, * 2 ' i*2’ ^S' ^g)» an< * 

bi(z) is a function of the unknowns Pi j , i=0,l,...Ng 

j=0,l, . . .F-i-1 

that we need to determine in order to solve the linear system (3.2.2) 
for the Pi(z). Using Eqs. (A.l) through (A. 3) of Appendix A, the 
Pij may be expressed in terms of Poj, j»0,...F. To find these 
Poj , we need (F+l) equations, obtained as follows: 

. Equation (A. 4) solved recursively yields Ng equations 
. One equation is obtained by Pij » 1 

. Considering the Cramer solution Pi(z) » ^t^izj where 
Ai(z' is the matrix A with i-th Column replaced by vector b(z). 
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It can be proved that det A(z) has Ng unique roots 
2 ^ in (0,1). The existence of the Pi(z) requires 

det Ai(Zj) ■ 0 for j»l ,2 , . . .Ng (3.2.3) 

which leads to a set of Ng equations. 

At the costof great computational effort in this intermediary step, 
one can manage to solve for Pi(z) . Related details can be found in 

FIS 7? ! and related work by Chang. 


a) PS traffic performace 


Case B finite: In addition to the average packet time delay, 

an appropriate measure, useful for buffer size design, is the 
probability of overflow. 


P o “ Pr £i+j « B+FJ - EjJq Pi, B+F-i (3.2.4) 

Case B infinite: Pkt performance is measured by the average 

time delay 


E(T) - i 


: j s 

L i»j 


d Pi (z) 
dz 


(3.2.5) 


Only simple examples (Ng*l, Ng“0) have been treated in the 
literature [wei 80* using this method because realistic cases 
involve unmanageable computations inherently encountered in 

recursive solutions and finding the roots of det Aft) = 0 in 


the disk (0,1). 


b) CS traffic performance 

The probability of steady state occupance of the M/M/Ng/Ng 
system governing the CS process is 
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, where p^ 


(3.2.6) 


The blocking probability, the CS performance measure, is then 
a special case: 


Pg * ^Njj slots are busy^ 


N< 

p l 

"TT 


/N s ! 


(3.2.7) 


m , 


c i /r!! 


It should be noted that the above expression for Pg, known in 
telephony as the Erland-B formula, is a approximation in this 
case because CS customers access the queuing system only at the 
start of the frame. Therefore, the Erlang-B formula is valid only 
for short frame durations. Fig. 3.4 shows the dependence on 
P D and the validity of the Erlang-B formula with frame duration. 

D 

In conclusion, the model, introduced here for completeness, 
highlights the computational complexity, mentioned earlier that 
is involved in analyzing integrated switching. Consequently, 
this model is abandoned for the rest of this report. Instead, 
approximate methods to analyze PS traffic performance are 
considered in the next section. 


3.2.2 Approximate analytic models: 

a) Formulation as M/D/1 queue with ON-OFF server 

In the frame configuration that we have been describing 
so far, we seperate the CS region (Ng slots) from the PS region 
Ng slots). However, as far as implementation is concerned, the CS 
slots are not confined to occupy only the CS portion of the frame, 
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but can be spread out or interleaved along the frame in a deterministic 
or quasi-random manner. Thus, a packet arriving at any instant within 
the frame may attempt service at the start of the next slot if this 
latter is available (not occupied by a circuit). So, a time slot 
selected randomly in the frame is intermittently available for 
packet transmission. Clea^’v, the probability of the slot availability 
at some time, denoted by c depends on the frame management policy 
(fixed or moving boundary) and the num ier of circuits in progress 
at that time. Kith this frame implementation, the system lends 
itself to be modelled as in Fig. 5.5 by an M/D/1 queue where the 
server is available for PS customers with probability o given by: 


* Fixed boundary 



V n b 


(3.2.8) 


* Movable boundary 
a = 


F-e 


(3.2.9) 


where e is the average number of slots occupied by CS customers 
within a frame. 

An expression for the queue length of the randomized service 
M/D/1 system may be obtained as in[0CC 77] . 


Let Pn = Pr in pkts are present in the system at the start 
[of the current data slot. 

Y k “ Pr new arrivals during a data slotj 
Then , as in joCC 7?J , 


Pn 


E 


1-cOP 


n-k 


+ o P 


n+l-k 


k_ 


'n 


'V aP ! : 


(3.2.10) 
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and the moment generating function is 


P(z) 


aP 0 (z-l) 



(3.2.11) 


. P Q * Pr ^empty systenij 

* 1-X,A* = 1-X, — (3.2.12) 

L 2 a 

t* is the effective constant 
service time (slot duration) 

. f!(z) = '<;£* is the mgf of the pkt arrival (3.2.13) 

• u J process. 


Using this expression, the average number of pkts in the integrated 
switching system is then 


l- 


z = 1 


(3.2.14) 


Assuming that packets arrive according to a Poisson process, 
and using Little's result, we get the average packet time delay 

A (l-P 2 /2) 


L_ 

■ i 


O-D- 


(3.2.15) 


or normalized 

I 

A 


1-P 2 /2 
a - oZ 


(3.2.16) 


It can be noticed that this result is similar to the one of a 

A(l-P 2 /2) 

normal M/D/1 queue (given by — *— ). Clearly, the effect of 

i-p 2 

service interruptions, which occur with relative frequency o, 


increases the delay suffered by an average oacket and limits the PS 
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traffic Intensity p 2 to be smaller than a, otherwise the system 
gets unstable. 

This model Is used In the next chapter to evaluate PS traffic 
performance in an integrated multiplexor transmitting over a satellite 
channel . 

b) Formulation as a single server subject to "Breakdowns" 

Contrary to the previous implementation model, the CS slots are now 
assigned adjacently in the circuit region of the frame (N<- slots). Thus, as far as 
the transmission of packets is concerned, the FFFB or FFMB frame structures can 

be modelled as a single -service system whose service is interrupted during the 
portion of the frame occupied by CS traffic (see Fig. 3.6). A packet arriving 
during the PS portion of the frame starts service at the start of the next slot, 
whereas a packet arriving during the CS portion of the frame "sees" the server 
in "breakdown" state and thus must wait until the start of the PS portion to 
be served. (See Fig. 3.6). Packets are assumed of constant length. 

As before, the CS slot occupancy is characterized by the probability 


and 


= Pr k occupied CS slots 


P*/k! 

y N'S m 

^m-O p 1 


/m! 


(3.2.17) 


Pr blocking 



(3.2.18) 


Here, o, as defined previously, will represent the fraction of the 
frame reserved for PS traffic (and 1-a that reserved for CS traffic) . 

We focus on PS traffic performance analysis by evaluating the 
mean time delay using a "work" concept and taking into account the 
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"repair time of service breakdowns". [Vu 8lj . A typical realization of 
the "work" function is portrayed in Fig. 3.7. 

Let t = AF : duration of the frame 

A * pkt service time (or slot duration) 

T » mean system time of a pkt from arrival until departure. 

A new pkt may arrive during the CS or PS portion of the frame with 
probability (l-o) or o, respectively. 


Case of fixed boundary 

Assume first that a packet arrives during the PS portion of the 
frame, and let be its system time. The overall time spent by 
a packet in the system consists of three components, namely: 


. ~2~ : expected time before pkt service will be interrupted. 

. ( 1+ x 2 T i ) A : Residual work to be completed on average by the 
packet before leaving the system. The term 
*2^1 re P resents the average number of packets 

that would arrive during the time (T.) the 
packet is in the system. This is equivalent 
to the number of packets that the arriving 
customer will find on average in the queue, 
and which have to be served before the considered 
packet can go to service. 

"" (1 + A,T.)A - C t/2 ~ 

. t—1 (l-o)r: Service repair time which is 

0T equal to the number of breakdowns 

times the average interruption 
period. 

where x^ , x denote the smallest integer > x and the largest 
integer < x respectively. We can approximate x by x+1 | and 
therefore the mean system time 6f a pkt arriving during the PS 
portion of the frame is 
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*1 • r- ♦ ci** 2 t,)4 *[- 


(l+X^^A+o t/2 


OT 


J 


(l-o)x 


(3.2.19) 


If the packet arrives during the CS portion of the frame, it 

will be delayed, similarily, by 

. (l-o)i : expected time until start of PS portion (start 

" of service) 

. (1+X,T 7 )A : residual work (T, is Average System time) 

w w te 

(1 + A ? T->)A 


U + a 2 i 2 ja 

OT 


(l-o)x : service repair time 


therefore, we have 

. T 


1+A,T,)A , 

T ? = (1*0)^ + (1 + A 2 T 2 )A + — ■■ 


(1-o)t 


(5.2.2U) 


The total pkt system time is the weighted sum of and T-,, 

T - + (!-o)T 2 


(3.2.21) 


If the number of service interruptions (breakdowns) is large, 

we can approximate J_xJ = x in Eqs. (3.2.19) and (3.2.20). 

Thus, Eq. (3.2.21) becomes 

, | (1+X,T,)+Ai I |(l+X ? T,)Aj 

T * (1 -oy\ ♦ (l+XjT^A + ^ Z -J (1-o)t + |— g T ^ JCl-a) 

(3.2.22) 

Using Eqs. (3.2.19) and (3.2.20), we get the normalized system 
time [yu 8l] . 

T d-o) t/ 2 • i 


1 - X 2 4 


(3.2.23) 


Case of movable boundary 


The effective boundary is represented by the normalized average 
number of busy slots in the CS region, that is 
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N- kP t 

E(Y) " E k«o » where p k 


and its second moment is 


7 *S 1 P k 
E(Y " ) “ Z k=0 k “ 77 
N S 


p k /kl 

“Hi — 

Co p T /m! 


(3.2.25) 


(3.2.26) 


The average system time delay encountered by' a packet is 
estimated in the same manner as in the previous case, noting that a 
packet that arrives during the CS or PS portion of the frame is 
expected to wait E(Y*')T/2E(Y)or Efl-Y)" t/2 E(1-y) before the CS 
or PS portion ends, respectively YU fclj . 

The partial system times, conditioned on the arrival time 
within CS or PS portion, are respectively 


cr 2, [ (1+A,T,)A| 

W * (1 * X 2 T 2> i * hC^-J E(Y)I (! - 2 - 27> 

(UA-TJA- E(l-X) 2 t/2E(1* ') -| 

L 1 E(y)t 

(3.2.28) 


E (i- Y) 2 + (1 + X,T,)A + 

2E(1-y) 


E(1-y)t 


The total time is 

T - [l-E(Y)] T 1 * E(y) T 2 

After some calculation, we get[YU 8lj. 

at . 1 

I . JE T TT7T 

A X 2 A 

1 ’ nm 


where 


oi * 


E(Y 2 )^E(Y)-2E(|)EjY 2 )-2E 2 (Y) +2E 3 (y) 
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3.3 INTEGRATED SWITCHING IN A MULTIPLE BEAM SATB.LITE SYSTEM 

In the o on text of the NASA 30/20 GHz studies, the problea of multiple 
eooess and resource allocation has been previously considered in a packet 
switched traffic environment [STE80] . The same problem is adressed here in an 
integrated circuit-switched (GS) and paoket-switohed (PS) traffio environment 
arising in a heterogeneous population of users to be aooommodated by a 
multibeam satellite system. We first conoern ourselves with the integrated 
station access problem. Several schemes are then presented and their 
performance compared. 

3.3.1 Uplink Capacity Allocation 

The satellite system is assumed to serve N zones via uplinks and 

z 

downlinks of the same capacity (e.g., 2 Mbps). The same average traffic load 
flows between zones. M earth stations located in each zone collect traffio 
from stream (CS) and bursty (PS) terminals in that zone. The eooess problem 
of these terminals to their earth station is not considered here. Therefore, 
earth stations within a zone are viewed as integrated users, sharing the same 
uplink channel in a TDMA mode. The oapaoity (C bps) of the satellite channel 
is represented on the time axis by a fixed frame structure managed on a fixed- 
frame, fixed-boundary (FFFB) or fixed -frame, moveable-boundary (FFMB) polioy. 
The frame is divided into F time slots of duration A . A nominal boundary 
separates the CS portion (N s slots) from the PS portion (N 0 slots). The 
multiple access of stations within a zone as illustrated in Fig. 3*8 leads to 
the following resource allocation problem. 
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Capacity allocation at the station aooaas level: 

Qivan a aatallita uplink ohannal of oapaolty C bps, abarad by N CS/PS 
integrated users in a TDHA node, we seek strategies to assign this oapaolty to 
the integrated users in order to aohieve the following objectives: 

e Efficient utilization of the bandwidth (maximum utlisation of the 
frame). 

e A desirable level of blocking for CS traffio (e.g., Voioe Calls), 
e The best possible delay-throughput charaoteristio for PS traffio. 

(l.e. , minimize delay, maximize throughput). 


3. 3. 1.1 Circuit Capacity Assignment 

Focusing first on the CS capacity (N slots) assignment problem, several 

s 

schemes are now presented (see Fig. 3.9). 
a) Fixed assignment. 

N 

The same number of slots, N « is assigned to each station 

during eaoh frame in a fixed manner as in Pig. 3.9a. The olrouit arrival 
process from each station is assumed Poissonian with the same rate X and 
the cirouit holding time is exponential with mean Assuming that the frame 

duration is short enough ( compared to circuit holding time), the aeoes of CS 
customers from a given station to the channel is governed by an M/M/N/N system. 
Again, the state probability of this system is 


P 


k 


o |Vk! 

m»0 


(3.3.D 


where 
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These CS customers arc b looked whan all N slots art busy, that la with 
probability. 


pf/kl 

P r lk .Ml . P B . — ji- (3.3.2) 

V _®/_f 

o-O P l' 

A quantity or crucial interest, used later in paoket performance 
evaluation, is the total average number of slots oocupled by the CS oustomers. 
It is given by 


M 

e * M z k P, 
k-0 * 

z M Pl (]-P B ), (3.3.3) 

after some calculation. 

It is clear that this fixed assignment scheme leads to bandwidth waste during 
low traffic oeriods because many slots remain unoccupied. 

b) Shared-pool assignment 

The corresponding frame configuration is illustrated in Fig. 3.9b. 
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A pool of U unaselgned slots Is sods svallabls In tbs CS portion. Tbs 
remaining (N g - u) slots srs ssslgnod to tbs M ussrs In n flxsd manner. Ths 
froo pool U is sbsrod by all CS oustomere on a de man d basis, that is if ths 
N slots ars all busy than ths overflowing traffio has aooass to ths pool. 

The probability of overflow is just 


V 


pJ /n : 


? ■/ . 

»£o °i / b; 


(3.3.4) 


The intensity of the total overflowing traffio competing for the pool is given 
by 

p u a M °1 P 3 1 (3.3.5) 

The overflow arrival process is no longer Poisson, but under heavy traffic the 
Poisson property remains approximately valid. With this assumption, an 
overflow CS customer which finds the pool full is blooked with probability 


P B ■ 


(P u )“/u.* 


! 


(3.3.6) 


m-0 <°u 


Finally, the total blocking probability for a particular CS customer is given 
approximately by the produot 

p„ . r t , h (3.3.T) 

The average number of busy slots oan be derived as follows. 


» « p i<i - v ♦ 1 - y 


c * M P x v1 - Pg) 


(3.3.8) 


o) Demand assignment. 

In the scheme of Pig. 3. 9c . any arriving CS oustomer is assigned a time 
slot (if available) in the CS portion of the frame, on a demand basis. The 
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intensity of arrivals is Mp^ and P B is, similarity, glrcn by 
N N # 

P B . (M Pl ) */N s I / l M Pl ) B /«l (3.3.9) 

m»0 

The average slot ocoupany is again 

e * - P B ). (3.3.10) 

d) Shared Channels Assignment 

In the previous model, the pool of U slots is of arbitrary size and its 
access is unrestricted. Beoause of traffio variations, a particular station 
may occupy most of the pool, creating congestion among other ineoming CS 
channels. To avoid this, the same model is refined in the following way: The 
pool is divided into R shared channels consisting of slots each. Any 
station, on demand, may take one of the R channels (if any) and use up to 
slots. The station must release this channel when it is not needed to make it 
available for the other stations. The corresponding frame configuration is 
illustrated in Fig. 3.9<t. It should be noticed, as opposed to the previous 
model, that N,R,N A are variable parameters and are related by 

N s * MN ♦ RN a (3.3.11) 

It is of interest to find optimal values of R, and N (or equivalently, 
optimal pool size u * RN^) which yield minimum blooking probability for an 
arriving CS customer. 

The queuing system handling the M CS stations with their assooiated 
shared channels can be viewed as a set of interoonneoted M/M/m/m queues (where 
m is either N or N ♦ N^). In Appendix B, we analyze this system to find the 
state equation from which we derive the blooking probabilty Pg for CS arrivals 
as well as the average state oocupancy. 
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Choice of R end N^: 

By examining the expression for P B , one would like to ohoose values for R 

and N (or equivalently N^), say R* and N* (N^ # ), that reduoe the blocking 

probability as much as possible. To be able to choose these values, one needs 

to know how Pg varies with R and N(N^). The search for (R*,N # ) was carried 

out numerically for several examples. It was observed from these examples that 

N s 

P. decreases when N increases for R fixed (0 < N < — -) , and that the smallest 
Pg is attained for N # * ■ M a - -1. On the other hand, among the minima of Pg, 
the smallest value tends to correspond to the value of R given empirically by 



Prom the equation 
empirical results 


N g s N*M ♦ R*N^ we deduce = 2 using the general 

M 

N* s — 1 and R* s — stated above. 


(e) Impact of digitization rate on circuit capacity assignment. 

It was mentioned earlier that the services to be integrated in the 30/20 
GHz system occupy a large spectrum, from narrow to wideband signals. The CS 
portion of the integrated multiplexor TDMA frame must thus handle transmissions 
with a variety of digitization rates (bits/sec). This suggests the idea of 
assigning channel capacity (time slots) in proportion to these rates. Let us 
assume that, classified by decreasing order of rate, we have K clases of CS 
customers. Frame slots may be assgined as in Table 3, where 


[ *1 > R 2 “• > \ 


°2 > “k 


(3.3.12) 


and 


a - 


R i 


SR. 


(3.3.13) 
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Table 3.2 Traffic Classes by Rate 


Class 

Traffic 

Intensity 

Rate 

Slot Allocation/frame 

1 

O 

II 

*1 

slots (bytes) 

2 

C 2 

*2 

a " 

2 

1 

t 

» 

t 

1 

f 

I 

t 

1 

t 

t 

f 

« 

• 

» 

♦ 

i 

1< 


R, 

i «„ 


k 

k 

k 

J 


Assuming that the frame duration is much shorter than the average circuit 
holding time (exponential), this system can be viewed as a K-dimensional 
Markov Chain. The state equation is given by the closed form expression. 

p"i 

P(n 1 ,n 2 ,...n |c ) = P(0) n “n^T (3-3.14) 

The constant P(0) is determined by the normalization condition 

Z P(n) = 1 (3.3.15) 

S 

where S is the space of all allowable states, defined by the resource 
constraint 

T 

2 H = ♦ c 2 n 2 ♦ 

From Eq. (3-3.15), we have 


• • • 


♦ Ws 


( 3 . 3 . 16 ) 
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P~ l (0) 


"1 

E 


2 

E 


n l"° n 2*° 



k 

n 

i-1 
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(3.3.17) 


where is the maximum number of class 1 customers In the system, Ig in the 
maximum number of olass 2 customer given that class 1 customers are in the 
system, etc. That is, 



where [xj denotes the largest integer <_ x, 

A customer of class i is blocked with probability 


(3.3.18) 


P n * r P(n) (3.3.19) 

B i !> 

where the sum is over states such that n e S and n ♦ 1 S, where e^ is the 

i th call entry vector (all zeros except 1 in the i^ 1 position). 
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Exponential Assignment: 

Assume, without loss of generality that 

( 3 . 3 . 20 ) 

where a is an integer constant greater than 1, and m is an integer. Time 
slots are assigned to a class i arriving customer according to the rule 


M s « a*" 1 * 


k-i 

a ^ = a slots 


In this case, we have 


( 3 - 3 . 21 ) 


P" 1 ( 0 ) = z 


l ... z 


k 

n 


n ^-0 0^*0 


n -0 i*l 
k 


n, 

°i 1 

n i : 


( 3 . 3 . 22 ) 


where 






v k“l 

N s j* 

a l S a "* 1 

k-1 k-1 

a m - a n. 


= m 


k-2 


c^m-n^ = a(Ij-n^) 






( 3 - 3 . 23 ) 


k-1 k-1 . k-2 . k-i+1 . 

a m - (a n^+ a 0^+. ..+a n i-l^ 

k-I 


* °[“ 1 ' 2 * ai '\ ♦ — *"!.!>] * “t'l-rVi] 
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This problsa Is known in olassio ts lap hone trunking as the multi-ohannel 
problem beoausa of tha N s slots (or channel) wittun a fraaa. It arises 
typically in a T1 transmission system in which, for example, 64 Kps and 32 
Kbps rates are to be handled. 

For illustration, let us consider the case of K s 4, a = 2, N * 2^m. 

The N slots are allocated to arriving customers (calls) as in Table 3 *3. 
s 


Table 3.3 Slot Allocation 
k = 4, a s 2, N s = 8m 


Class 

Slots allocated 

1 

2 3 s 8 

2 


3 

2 1 s 2 

4 

rH 

II 

O 

CM 


We are interested in the blocking probability for each traffic class. For 
instance, this probability can be computed for class 3 in the following way: 

P B ‘ 1 } 1 


£p(n^,n2,n2,8m-8nj-4n2-2n2)+P(nj,n2,n2,8m-8n^-4n2“2n2-l)J 


+ P(0,0,4m,0) 


(3.3.24) 
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The total b looking probability la given by the weighted sum 


k-4 

P. * E 


B 


1-1 


k 
r p 
J- i 


(3.3.25) 


j 


The average number of busy slots (channels) in the system is given by 
E(8nj+4n 2 ♦ 2n 3 ♦ n 4 ) * 8E(n l ) ♦ 4E(n 2 ) ♦ 2E(n 3 ) ♦ E(n„). 


And it is known from before that: 


E( n i> 8 p i (1 " (3.3.26) 

i « 1, .... K = 4 

which is Just the average number of class i customers in the system. 

Consider the simple numerical example of K = 2, a = 2, N g = 32 = 2 x 16 
= 2 m. From the previous expressions we have 

-1 

“1 

°r °2 2 

P(n,.iu)s* r E -fr ~ tt I r ~r zrr~ (3.3.27) 


C m 2(m-i) pj p j 2 ] p"l P^2 

:n.,n ) = ( E l — 7 — ~Ti I Z~r — ~ 

1 2 ^ i=0 jat0 i. J. J V n 2 ’ 

= E jV(n lt 2m-2n^-l) ♦ P(n lt 2m-2n 1 )^ + ?(m f 0) 


m 

P- = E P(n., 2m - 2n.) 
B 2 n x -0 1 


03.3.28) 


(3.3.29) 


E(2n 1+ n 2 ) = 2 p 1 (1-PB1 ) ♦ P2 (1_P B2 ) 


(3.3.30) 


Considering the traffic mix 


8 * 



(3.3.31) 


as a parameter, the blocking probabilities for classes 1 and 2 are plotted 
versus the offered load in Figs. 3.1° and 3.11 for two values of B, 0.5 and 
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0.75. Similar paramatrie eurvea oan ba obtainad for k > 2 olasaes of 
oirouit-awitohad sarvioes with various ratas. Thaaa ourvas may form an 
important tool for tha ayatam daaignar so as to avaluata a dasirad dagraa of 
aarvioa (hara P 0 ) under a given traffio mix ( 3 ) oondition. In a sense, 
ranking tha users' traffio according to tha ratas is analogous to assig ning 
partial priorities. Figs. 3.10 and 3.11 show tha gap in the degree of service 
(blocking) between the two classes of users. This gap, as expected, is 
relatively more accentuated for light traffic. 
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3. 3. 1.2 Packet Capacity Allocation 

The fra aw any be shared by PS traffic under both fixed and moving 
boundary polioies. The average time delay, the paoket performance measure, is 
evaluated using the approximate rode Is described in seotion 3*2.2, namely: 
•M/D/1 with random server 
•Single server subject to "breakdowns.” 

The expressions for paoket time delay are given by Eqs. (3*2. 16), (3.2.23) and 
(3*2.30) Because of its superiority in performance, the FFMB scheme is given 
more emphasis here and is considered in the following numerical example. It 
is important to stress the fact that the two models above correspond to two 
different frame implementations. Consequently, these two models cannot be 
compared on a numerical basis. In the example treated in the next section, we 
used only the first model to evaluate packet performance . Clearly, the 
evaluation of this performance is straightforward using the second model by 
referring to Eqs. (3.2.23) and (3.2.30). 
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3.3.2 Performance Comparison 

The uplink satellite ohannel (assumed of capacity 2 tops) is 

synohrcnoualy elooked into time slots of duration A msec, grouped in frames 

of F slots. Data packets, assumed of oonstant length P * 1000 bits, are 

transmitted one per slot, so that A * - » ^9. * 5 msec. The oirouit and 

2 x10 

packet performances are now evaluated on the following example: P * 50 slots, 
N a * 24 slots, N 0 « 26 slots. 

Fla. 3.12 shows the variation of the blocking probabilities with the CS 
offered load per user, , for each of the four slot allocation schemes. We 
note, as expected, that the performance of the fixed assignment scheme is much 
worse than the two demand assigned schemes for small blocking probabilities. 
The shared channel (where the pool is accessed two slots at a time) scheme is 
somewhat better than fixed assignment for this uniform traffic environment, 
but still considerably worse than the demand assigned schemes in the same 
range of traffic. 

The comparative performance of the schemes reverses in the case of 

packet delay, however, since in demand-assigned schemes, fewer slots remain 

unused. Fig. 3.13 showing packet delay as a function of paoket-switched 

traffic intensity c> 2 , with ^ • 2 Erlangs eireuit-switched traffie per user, 

indicates that the fixed assignment scheme is now the best, although not by 

N B 

very much. Note that the curves here are shown dashed beyond p 2 * — |r > 

0.52. This indicates that Eq.(3.2.l6) on which these curves are based, is 
presumably no longer valid beyond that point. This is because in the original 
model [0CC 77] , it was assumed that the CS and PS processes are uncoupled 
which is only valid approximately if one does not attempt to use the added 
slots made available by the movable boundary to drive the system with a PS 
traffie intensity e> 2 > 0.52. The added slots may be used 
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to reduce packet time delay, not to increase throughput. 

All four schemes utilizing a movable boundary strategy show considerable 
reduction in paoket time delay over a fixed boundary strategy. For this range 
of parameters, then, it is the movable boundary that provides the baslo 
improvement in paoket performanoe. The added effeot of the speoific circuit- 
switched assignment strategy represents a relatively minor perturbation. In 
Fig. 3>16 we observe that the performanoe of the demand assignment strategy 
deteriorates with increased circuit-switched load more sensitively than the 
other strategies. However, this effeot must be interpreted with care sinoe 
blocking probabilities become unacceptable for high range of c^. 

The comparative study carried out here may be extended to the case of a 
more heterogeneous-traffic user environment, where results may differ from 
those obtained in this work. 


3.4 ON BOARD PROCESSING AND STORAGE 


Many advantages are associated with equipping the satellite with pro- 
cessing repeaters Instead of conventional translating repeaters. The 

former provide on-board demodulation and remodulation which permits switching 
to be performed at baae band frequencies, resulting in leas hardware complexity. 

The effects of noise and Interference can be reduced by appropriate control 
and regeneration of the detected data. 

More importantly, the availability of detected data (bit streams, packets, 
messages...) at the satellite, with the addition of some intelligence and 
buffering capabilities, suggest the use of the satellite as an Integrated store 
and forward node performing the sequencing and switching of the CS/PS traffic 
from uplinks to downlinks (see Fig. 3.17). 

3.4.1 Integrated Downlink Channel 

The outgoing pipe consists of N downlink channels connecting independently 

z 

the satellite to the N zones as in Fig. 3.18. A buffer is associated with 

z 

each downlink channel and the collected traffic is routed to the appropriate 
buffer corresponding to its destination. The channel is assumed accessible to 
CS traffic with preemptive priority over the PS traffic. CS traffic may be 
managed on a blocking basis. 

In this section we study the buffer behavior and the resultant queuing 
delay at the satellite level. Ue first assume that both CS and PS traffic 
have the same statistics, namely Poisson arrivals and exponential service time. 

A frame structure is used for the entire satellite-earth pipe and each downlink 
channel is assigned a sub-frame of 2 slots , The first slot forms a CS subchannel 
and can be used, if free, for packet transmission. The two slots are not necessarily 
of the same duration. Thus, the Integrated channel can be modelled by a single- 
server queuing system in which the capacity of the server alternates between 


two levels In a Polasonlan manner. This la an application of a general queuing 
modal p recanted In Qyec 60 . 

Next, we Introduce priority traffic claaaae aaaodated with ranked services. 
Following the work by Shantlkumar JjSHA 8lJ , his model Is extended below to our 
application of an intagrated-switched channel, assuming constant packet length. 

3. 4. 1.1 Two-Dimensional M/M/1 System 

As far as PS traffic is concerned, the server (channel) oscillates between 
two states denoted 1 (channel occupied by a circuit call) and 0 (channel free). 
The time that the server remains at a given state is exponentially distributed 
just as interarrival and service times of the CS process have an exponential 
distribution. The queuing model is shown in Fig. 3.15 and the corresponding 


state diagram in Fig. 3.16. 


The joint steady-state probabilities (P^ , 1*0,1; j»0, 1,...} 

are related 

by the following balance equations: 


*10 - "2<»*11 + *1*00 

(3.4.1) 

Q 2 < 0 >^ *00 -"2 (0) *0l*"l*10 

(3.4.2) 

Q2«)-2<l)-0 'll ‘ U 2 < H*1 ,j+l + ¥‘>'l.j-l * *1*03 

(3.4.3) 

[> 2 (0).y 2 (0)«Q P o] - u 2 (0)P o>j+l ♦ X 2 (0)P o j .j + VlJ 

(3.4.4) 


In the case of finite buffer size (say B) , the system of equations above, 
together with the condition - 1, can be solved recursively for the 

P^s. This involves (?B+3) equations, which may be manageable if B Is not very 
large. Quantities of Interest, such as average buffer content or probability of 
overflow are readily derived when the P^'s are known. This case fits very well with 
the a Plication to the on-board buffer for a permissible probability of overflow. 
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If « large memory capacity le provided on the satellite so that queue 
lengths cannot exhaunt the storage space, we can aasume infinite buffer site. 
This case is treated In the following. 

Since the average persistence time of the server (channel) In state 
0 and 1 la 1/Xj and 1/y^ respectively, we have the probabilities: 


- Pr 

II - Pr 

0 


E 


server In state 


G 


serve*.* In state 



Jll 

1+r 


1 


1 

l+o 


1 


(3.4.5) 

(3.4.6) 


The average arrival and service rates for the packets in this queuing 
system are defined as follows: 


x 2 - n 0 x 2 (0) + Hj x 2 (i) (3.4.8) 

W 2 • n Q u 2 (0) + IIj u 2 (l) (3.4.8) 


It should be noticed that the packet arrival process Is totelly independent 

of the server state, i.e., X 2 (l) ■ X 2 (0) ■ X 2 , but the notation above 

Is kept for convenience. It will eppear later that a proportionality relation 

X 2 (l) X 2 (0> 

between the traffic parameters - ^ ^ ^ leads to extremely simple solvability 

of the problem. This relation can urn generalised and It Is discussed leter as 
a Flow Control mechanism. 

The set of difference equations (3.4.1) th* ough (3.4.4) can be reduced 
to the single equation 


x 2 (i)P lj + x 2 (0)P oj - 


u 2 (l)P l,j+l 


V 0)p o,j+i 


(3.4.10) 


Put in vector form* this becomes 
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with the vector definition 

Pj -[p^, P Qj ] and A 2 -[ x 2 (1) ‘ Y 0) ] (3.4.12) 


Eq. (3.4.10) recalls the M/M/1 result XP^ • yPj + ^. Thus, the queuing 
system under consideration can be viewed as a 2-dlmensional M/M/1 system. The 
average residual service capacity obtained from (3.4.10) by summing over j, 
is given by 

V 2 - T 2 - P 10 U 2 (1) + P 00 li 2 (0) (3.4.13) 

For the system to be stable, this residual capacity must be positive, i.e.. 


Recall that ^ and U 2 are given by (3.4.8) and (3.4.9). Now, defining the 
conditional m.g.f. of the system as 


G i (z) - Pjj 2 ^ , i - 0,1 
We can transform Eqs. (3.4.1) through (3.4.4) into 


(3.4.14) 

the following linear 


system in G^(z) and Gq(z): 


Aj(z) 

- X 1 


Gj(z) 


B j(z) 

" U 1 

V z) 


G o (z) 

- 

B 0 (.) 


where 


Aj(=) - j 

^ 2 (1) + Uj + W 2 (l7| 

- 

Q 2 (1)z + u 2 (l)/zj 

A 0 (z) * 

[X 2 (0) + Xj + u 2 (0^ 

- 

Q 2 (0)z + u 2 (0)/zj 


S 1 (z) •» u 2 (i) (~) P lo , for i « 1,0 


(3.4.15) 

(3.4.16) 

(3.4.17) 


(3.4.18) 
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To solve this system, we need to determine P^, i ■ 1,0 l.e. , the prob- 
abilities that the system is empty of packets in one state or the other. For 
this, we need two equation in P^q and Pqq- One equation is given by Eq. (3.4.13) 
and the second is obtained as follows: 

The Cramer solution of the system of (3.4.15) is given by 


G ± (z) 


det M i (z) 
det M(z) ’ 


(3.4.20) 


where M 


Aj(z) - Xj 

~ u i V z) 


and is M with the i th column replaced by vector B(z) 

It is known that det M(z) has a root at some z q , (zq| < 1. Therefore, 
det M^(z) must have the same root for a solution to exist. 


det M(Zq) * det M^Zq) * 0, for |z q | £ 1 


(3.4.21) 


This condition leads to finding the root of the third-degree 


polynomial QfEC 650 . 


g(z) - Az 3 + Bz 2 + Cz + D 


(3.4.22) 


where 


A = X 2 (l) X 2 (0) 


(3.4.23) 

b - - Q 1 x 2 (0)+x 1 x 2 (i)+x 2 (i)x 2 (0)+x 2 (i)u 2 (0)+\ 2 (0)ij ? (ir] 

(3.4.24) 

C = Qi 1 v 2 (O)+X 1 u 2 (l)+w 2 (O)y 2 (l)+X 2 (l)u 2 (O)+X 2 (O)u 2 (10 

(3.4.25) 


D * -u 2 (l)u 2 (0) 


(3.4.26) 


Combining equations (3.4.13) and (3.4.21), we can get the quantities 
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P 10 ,P 2 o|j reC wnich ® nabl ® us to obtain Gj(*) and G Q (*) from ( 3 . 4 . 15 ). 

The quantity of interest for us is the average number of packets in the 
buffer, given by 


“• ■ -3T [V 2 > + g 0 (2 3 

The final expression |*YEC 69j is 


E j»0 j (P lj + P 0j^ (3.4.27) 


z • 1 


p 2 r~ 

EL * u 2 (1)(u 2 (0)-X 2 (0)) P 10 +u 2 (0)(u 2 (1)-X 2 (1))P 00 

i- p. L 


- (u 2 (O-\ 2 (1))(M0)-> 


-y"D 


(Xl+Ui)(u 2 -X 2 ) 


(3.4.28) 


It is important to observe that the first term is the average queue length 

X 2 

of an M/M/1 queue with utilization factor p 2 * — — . In the case of an 

integrated channel, and during the long holding time of a CS call (— ^ 3 min.) 

u i 

the service capacity for packets is reduced (i^O) < y ? (0)) and the queue builds 
up to arbitrary large sizes. A sensitivity analysis would show that the second 


tern in (3.4.28) grows very rapidly with ratios such as 


u 0 ( 0) X, 


. The 


‘1 


‘1 


first ratio compares a call holding time to packet service time, the second 
ratio indicates the average number of packets arriving when the channel is 
(partially or totally) occupied by a circuit. 

This point is investigated in the following special case: For example, consider 

a single slot channel, i.e., u 2 (l) “ 0, with 

X 2 (l) * x 2 (o) - x 2 

y 2 (0) * u 2 


the parameters become 


x 2 * x 2 


1*2 * ” ^ 0^2 


(3.4.29) 
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From Eq. (1.10), w® have 


P B IT • ■ i - 

r 00 U 0 U, 


H 0 “ p 2 


(3.4.30) 


Substituting these parameters in Eq. (3.4.28), we get the average number of 
packets in the queue: 


EL 


7 2 


l - P« 


1 - P 


V 1 - V- 


where 


IT 0 u 2 


(3.4.31) 


(3.4.32) 


We readily observe that the ratio 


which is in the order of 10 to 


10** (e.g., typically = 100, y^ * 0.01 in Voice Data traffic), makes the 
second term much larger than the first term of Eq. (3.4.31). 


Using Little's result, the average packet queuing time is D 


EL 


. (3.4.33) 


This delay, which is the packet performance measure, is plotted versus the 
PS load in Fig. 3.22 for the same channel availability Hq but with two different 


ratios 


of 10 4 and 10 . The gap between the two curves is very significant 


and increases drastically with X 2 ' The packet buffer behavior in terms of 
time delay is depicted in Fig. 3.24 versus CS traffic load for a fixed PS traffic 
intensity (P 2 * 0.3). 


3.4. 1.2 Flow Control 

From the simple numerical case treated above, it appears from the curves 
of Figs. 3.22 and 3.23 that the queue of packets builds up during periods of high 
utilization of the channel by CS traffic during periods of peak circuit holding 
times, the accumulation of packets causes very long queuing delays if the stability 
condition is not violated. Thus, the need for a flow central (FC) and/or buffer 
management mechanism is apparent. This need is more urgent for our underlying 


application in tha satellite system because of long distance propagation mad on- 
board weight requirements. Next, a simple feedback control mechanism is presented 
to regulate Xj, the arrival rate of packets, according to the state of the 
channel (server). If a circuit is being transmitted (served), the feedback 
controller reduces >2 to X 2 (l), and if no circuit is on the channel the controller 
does nothing, i.e., X 2 “ X 2 (0) (see Fig. 3.19a). Interestingly enough , if this 
FC mechanism adjusts the packet rate X 2 such that 

X 2 (l) X,(0) 

— — rrr ” jnr. ■ 0 (constant) (3.4.33) 

then, it can be readily proved by induction QeC 60 that the general state 
of the system is described by the joint probabilities P j- j “ for 

i = 0,1. 

Therefore the packet queue length is governed by 
Pj * (Hq+II j) ( l-p)p*^ * (l-p)p^ which is nothing but the steady-state equation of 
an M/M/1 of factor p. It is very important to observe that the PS process is now 
completely uncorrelated from the CS process, i.e., the state of the channel 
(server) , and the number of packets in the system is simply 


where 


EL 


_e 

i-p * 


(3.4.34) 


*2^ m pkt arrival rate 
1*2 (0) Total channel capa. 


(3.4.35) 


The dynamic behavior of this oscillating queuing system is governed by 
the parameters characterizing the stastistics of the CS traffic emanating 

from users of various types. This traffic may be very dense (large X j) or 
very infrequent (small X^) with variable call duration or circuit holding time 
depending on the service (voice, file transfer, etc.). It is useful to qmantify 
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the behavior of the system under study when such extreme traffic situations 
occur. Two cases of particular interest, namely: 


• Very heavy CS traffic (X -*• ») and very short circuit holding time 


(Uj -*■ “) but with constant traffic intensity, i.e.. 


a (constant) . 


'1 


The system oscillates very rapidly between the two states and the 


probabilities of being in these states are 


n 


1 l+p. 


a 

1+a 


(3.4.36) 


n„ = 


0 1+c. 


1 

1+a 


(3.4.37) 


Equations (3.4.1) - (3.4.4) then yield 


P. . * P . , for i = 0,1,. 

13 oj’ J ’ ’ 


(3.4.38) 


From Eq. (3.4.13), we have 


P 00 *[^2 (0 > + 0 p i<°>) 


u 2 ~ X 2 


(3.4.39) 


p qo “ n 0 (1 " p 2 ) * (3.4.40) 

where __ 

— *2 

P 2 ~ ~ (3.4.41) 

v 2 

Similarily, 

P 10 3 ^l^ 1 ” ^2* (3.4.42) 

By inductive reasoning, we can show that 

P ij " V 1 “ p 2 } ^2’ 5 * 0.1,..., i - 0.1 (3.4.43) 

and the total state probability is 

Pj * (1 - P^p^. 


(3.4.44) 
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We observe thee the system converges to an M/M/l system with average 
parameters X 2 and u 2 . 

• Very light CS traffic (X^ -*■ 0) and very long circuit holding time 
(p^ -*•()), (e.g., transfer of a very long file in a single message) with 


The set of equations (3.4.1) through (3.4.4) becomes 


Xjd) p 10 - »j(D P u 

Q 2 (l)+u 2 (iT] Py - X 2 (l) Pj j.j ♦ Uj(i)P 1>J+1 


for i - 0,1 (3.4.45) 


which are the standard difference equation of an M/M/I system with 
X,(i) 

c i * p ~ (i) * Tn ot ^ er wor< * s » ^ X 2 (i) < y 2 (i) , i * 0,1, the system 
operates in two stable M/M/l modes with probabilities TIq and H^. 

Therefore, the total state probability is: 

P j * £ i*o,i n i Q ~ 0 2 (1)] oj(i) (3.4.46) 

Quantities of interest are readily derived from these probabilities. 


3. 4. 1.3 Generalization to Multislot Downlink Channels 


The control feedback mechanism described previously, which regulates the 
packet flow and uncouples the two traffic processes (CS and PS) , can be extended 
to the general case of a multislot channel where the frame is divided into 
F slots (Nj slots for circuit subchannel, N 2 for packet subchannel). The packet 
arrival rate is maintained proportional to the effective channel capacity 
P 2 (i) where 1 is the number of circuits in the system; that is 


X 2 (i) X 2 (i-1) 

u 2 (i) u 2 (i-l) 


* 2 (i) 


X 2 (i-1) 

F-(i-l) 


t * • • » 


N 


1 


F-i 


for i - 1 


(3.4.47) 
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X 2 (0) represents the uncontrolled pscket arrival rate (l.e. Max) and 
i^CO) the maxlmun capacity of the channel. 

We discussed this case in section 3.2 .1 . , where the packet subchannel 
was modelled as an M/M/S/® system with S » F - 1. Recalling the two- 
dimensional Markov formulation, it was not easy to find the state probabilities 
Pjj of the joint CS/PS stochastic process (i,j) and we saw how computationally 
complex the problem is. Fortunately, the FC scheme described above (see Fig. 
3.3.19b) transforms the problem to an extremely simple one because it uncorrelates 
the two processes (CS and PS) arising in the integrated channel, that is: 


P tJ * P t Pj (3.4.48) 

This fact can be shown by inductive reasoning on the set of balance equations 
Q'EC 60 just as it was shown previously for the case of Nj * 1, N 2 » 1. 
Furthermore, the marginal state probabilities are given by 

Pj/i! 

P i - — ^ — (II/M/Sj/Mj system) (3.4.49) 


. X 2 (0) X, 

V (1 - p)p 01/11/1 * th 0 -^( 0 ) * (N^H 2 )u 2 (3 - 4 * 50 > 


The performance of the system is evaluated by 

probability P„ and average packet queue length EL. 

N 1 

V /N l! 


E ni> 


the circuit blocking 
These are given by: 


(3.4.51) 



EL 


(3.4.52) 
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where 

X 2 

p " 0> 1 ’H> 2 >v 2 


(3.4.53) 


A simulation study would be needed to compare the performance of the system 
with and without the FC mechanism. It is mentioned in |kWO 8o] that a simulation 
experiment led to a performance improvement by a factor of 1000. 


3.4.2 Priority Strategies in On-Board Switching 

In the heterogenous environment of the satellite users that we described 
previous ly, an interesting issue is the assignment of priority classes among the 
data services. This priority ordering may be based on packet length, conversation 
duration, signalling messages, etc. We consider a set of K packet priority classes 
arriving at a given downlink buffer modelled as a single server (channel) queuing 
system. The availability of the server is randomly interrupted by CS traffic 
transmission. The system parameters are as follows: 


• Priorities 


# Poisson arrivals with rates Xj,...,X^ packets/time unit. 

• Packet service interval * one time slot ■ one time unit. 

# Pr ^erver available^ « a 

• FIFO service for same priority customers. 


Following the model in 


[SHA 8lJ , 


we define Che following quantities: 


0 waiting time: 
, service time: 


time from a customer's arrival to its first service 
attempt. 

time from first attempt to completion of service. 


. Queuing time: waiting time + service time 


The average queuing time EW f for priority-r customers is derived using level 
crossing analysis J^HA 8^ on the virtual waiting time V* process defined as the 
waiting time of a customer r if it were to arrive at t. If an r customer arrives 
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at t - 0 and tha buffer la found empty then It may attempt service; otherwise 
it has to welt until all higher priority customers who may have arrived ere served. 
Let G r be the random time needed to clear the buffer from these higher priority 
customers. So, just after t a 0, V p • and then drops llnearlly until it 
reaches zero or an r-customer arrives. In the first situation, an r-customer 
arrival Increases V* by G r because he either attempts service or has to wait for 
the higher priority customers to be cleared. In the second situation, the 
arriving r-customer waits until the preceedlng r customer Is served (call this 
service time) and the buffer is emptied of higher priority customers (call G^ 
this time). So V* jumps by the amount H r “ ® r + G^, ( see sample service realization 
and V* sample path [Iha 8l] in Fig. 3.20 and 3.21). 

When the buffer gets empty, a new cycle is initiated. We assume that the 
process V* converges to a steady state, i.e., 

lim t ^ V* - V r (3.4.54) 


and let V r (x) be the PDF of V , i.e., 

V r (x) « P r Q r < ^ (3.4.55) 

The level crossing analysis consists of evaluating the expected numbers of times 
the process crosses the level x upward (U(x)) and downward (D(x)), then 


E U(x) - E D(x) 


The analysis is carried out in 


Qha 87) 


and the details are omitted here. 


average queuing time in slots, for the r-prlorlty class is given by 


(3.4.56) 

The 


EW 


q(2-a) 

2 r o-9 r~b -9 7T 

_ lJ L_ r- i_ 


+ i 


x 

4-0 i 


(3.4.57) 


where 


(3.4.58) 
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and la in packata/alot. 

Using Llttla'a result, the total average number of packets In the buffer la 

EL - E W f (3.4.59) 

Some observations need to be made: 


• The first term In the expression EW^ Is the total waiting time 
(from arrival until successful attempt of service) which has the same 
form as the well known priority waiting time |~KLE 7 5*] . 


EW 

r 


E T. 


(i-e r )(i-e r _ 1 ) 


(3.4.60) 


with 





°1 


(3.4.61) 


Here E Tq is the residual service time given by 


E T 


0 


0 ( 2 = 0 ) 

2 


(3.4.62) 


in this case. 

• In the case of one traffic class (K ■ 1) , we have 


o + 2(1-X.) (1 - X./2) . 

py m — ■ — x — 

fcW l 2(a-X j) o - Xj 2 


(3.4.63) 


and interestingly, the first term is just the result derived by 
Occhiogrosso et al. in Q)CC 7TJ , using a completely different analytic 
tool. However, the addition of the term j makes this result more 
accurate because it takes into account the fact that a packet arriving 
within a slot can attempt service only at the start of the next slot. 

• More Importantly, EL, the total average number of packets in the 
buffer, is the same under priority service and FIFO service. This is 
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readily shown by computing 

EL - eJ b1 Xj EM t (3.4.64) 

which is 

(i-0'^(l-9 i -i') + ^J “ X Q2(,-X) ^ (3.4.65) 

where 

X - rj -0 X j (3.4.66) 

Now, considering a FIFO service discipline, i.e., one class of 
traffic with Intensity Aj ■ X, we see that Eq. (3.4.63) is exactly 
the same as (3.4.65) after using Little's result. 

Thij result is quite important if buffer space is a major constraint 
because we are accommodating priority services with the same buffering 


EL - Z \ ± 


resources. 
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. APPENDIX A 

TWO-DIMENSIONAL MARKOV CHAIN MODEL OF THE INTEGRATED TDMA FRAME 

The state translton diagram for the joint CS/PS process arising In an 
Integrated TDMA frame Fig. A. 1 Is given in Fig. A. 2. Taking into account boundary 
conditions, balance equations can be written by equating the rate of entering 
a state and the rate of leaving that state. This leads to the following set 
of steady stats difference equations £kW0 80]. 


• For i < Ng, we have 


[(X 1 +X 2 +iu 1 +ju 2 ^ p 1 j- j + < i+1 )M 1 P i+l j+X 2 p 1+l ^+X 2 p 1 ^_ 1 -»-(j+i) U2 P i (A.lj 


for j *0,1,... F-i-1 




for F-i < j < B + F - i - 1 (A. 2) 


[ V 1 V <F-I )ii2 ] f ° r 1 m B+F - 1 <*- 3 > 


• For i * N c , we have 


<VN S V^ 2 > P N s y + a 2 P M s .j-i- <J+1>ll 2 P »s, 


J+l 


(A. 4) 


for j - 0,1, . . . N -1 

o 




j+i 


(A. 5) 


V 1 


for j « 0,1, ... , 


A. 2 - 
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<N u,+N a y 2 )P M - x p >X P +X p 
S 1 B 2 N g J 1 N s -l,j 2 Ng^j_ 1 1 Ng_ l j+1 


(A. 6) 


for j * B+Ng 

Together with the normalization condition £ P 


ij 


1, Eqs.(A.l) through 


(A. 6) can be solved recursively to determine the steady-state joint probabilities 

P^. . However, for a realistic case, the number of equations Involved in the 

computations becomes quite large. This number is computed in the following way: 

For Q ■ 1, we have j ■ 0,1,..., B+F-i. This yields B+F-i+1 equations. 

Therefore the total number of equations is 

N s N +1 

Z (B+F+l-i) - ( — | — )(2B + 2N r + N. + 2) 
i-0 2 B S 

For illustration, an example with N g = 10, 5, B ■ 10 requires 231 equations 

to compute the , from which the quantities of interest in the performance 

criteria are determined. 


EFFECTIVE NOMINAL 

l BOUNDARY j BOUNDARY 
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STATE TRANSITION DIAGRAM 01 THE JOINT CS/PS PROCESS 
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APPENDIX B 

ANALYSIS OF THE SHARED CHANNELS MODEL FOR THE CS SLOT ALLOCATION 
PROBLEM 

The state of the queuing systaa presented in Section 3.1.1-d and 
depicted in Fig. 3.2(d) is described by the following closed fora expression 

M p"i 

P(n^,n 2 » . . . ,njj)-P(n) * P(o) H “ r > (B.l) 

X 1 

where 0, - and n. is the # of circuits from staiton i in the system. 

1 h 

The existence of P(n) can be proved by writing the following balance equation: 


(M * +M d u 1 )P(n 1 , ....n M ) = ^ Pfaj-l.i^, . . . ,n M > +P(n 1 ,n 2 ~l, . . . ,n fl )+ j 


M, nonzero terms 
a 


+ y ^ ,0^ »*•* *n^j) 4 , P(o^ ,02^*1 »••• * (B. 2 ) 

V 


^ 

H nonzero terras, 
u 


It can readily be shown by substitution that Eq.(B.l) satisfies Eq.(B.2). 
Now, it remains to determine the normalization constant P(o). This is done 
by summing P(n) over all possible states, and setting the sum to 1, that is: 

Z P(n) + E P(n) + . . . + Z P(n) - 1 (B.3) 

S 0 S 1 S R 

where state Sg means that K shared channels (ovt of R) are occupied. 
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The probability of the systes being in state S^, using synaetry, is 


Mi-M 

RTHi 


IF w ■ «i 

E P(n> - (") E 

Sk 


Oj-Hfl 


H+N. H 
E * E 

v "* 1 w° 


n *0 
m 


M P?i 

n P(o) -±-r 

i-1 V 


- P(o)(”) 



l * + *A p“i \ M / N 

1 -T. 1 11 2 

n^-N+1 “i* } i«K+l\n«0 





where 


. .M. _k M-k 

P(o)( k ) P q 


N+N, 


P n i 


n.-N+l n i ! 
1 


and 


n i*0 


±_ 

n. ! 


(B.4) 


(B.5) 


(B.6) 


Therefore, using Eqs. (B. 3) and (B.4), we get the following expression for 
p(o): 


P(o) 


R 

E 

k-0 




(B.7) 


Nov we proceed to calculate the probability that an arriving CS customer 
from one of the M stations (users) is blocked* This occurs in one of the 
two situations: 

• The station has all its N dedicated slots busy and all the R channels 
are occupied by the (M-l) other stations* 

• The station is using one shared channel (N^ slots) and all (N+N A ) 


slots are busy. 


-ab- 
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Without loss of generality, • custoaer f roa station 1 is blocked with probability 
Pg, obtained by suaaing over all blocking states. 


P - ( M_l ) 
B ( 0 ’ 


N 

E 

n 2 *0 


N 

r 


n M*° 


P(N+N^,n2 » • * • *n^) 




N+N, 


N 

E 


n 2 “N+l n^*0 


N 

E P(N+N.,.) 

V° 


H-l. 

+ <r-i* 


N+N, 


n ? =N+l 


N+N 
V A 

v ®* 1 


n 

z 

n R+r° 


N 

E P(N+N a ,.) 

v° 




N+N. 


n 2 -N+l 


N+N. 


n R+2 -N+1 n R+3"° 


^ P(N t » • * * >Rj||) 


n *0 

n 


(B*8) 


This can be written in a condensed form; 


where 


v * ("-*) 

B k«0 k 


k + - 4 


n+n a 

N+N 

N 

N 

E 

... E A 

V 

E 

n 2 *N+l 

n k+l -N+1 

n k+2*° 

n M*° 

1 

if k J* R 



0 

if k * R 




E P(N+k N A »n 2 , . . . *njj) 


1 

) . (B.9) 

d 


(B.10) 


and 


P(N+k N a , n 2 * • • • » n M )*P(o) 


N+k + N . 

fi t- 

<N+kV)! 

A 


n 2* 


M 


V' 


(B.ll) 


Using Eq.(B.7) for P(o), we get the final expression for P : 
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where 


and 


P m 

* cftH 1 ?) pV' 1 -" 

k*0 


(B.12) 

B 

T /^\ M-ffi 
E („)P , 

m«0 




^0^ +N A/(N+N a )! if 

k # R 


C(k) - 1 

. N 

^/N! if 

k - R 

(B.13) 





N+N c” 

p = l A — = — 

m = N+l m! 


(B. 14) 

N 





m-0 


m! 


(B. 15) 


The average number of slots occupied by customers from station 1 is 
calculated as follows: 


N 

E(n_)« E n. 
1 n^O 1 


R 

y .M-l, k M-l-k % 
E ( v )p q P(o) 

k-0 


V — 7 n I 

V 


N+N. 

E A n 
Rj«N+1 


R-l 


M-l x k M-l-k 


( i, ) p q 

k-0 K 


P(o) -±J 
n l 


(B.16) 


E(n ) R u i . 

1 _ r /M-1 n k M-l-k ' r n, 

1 ( fc )p q ; _ L _« 1 


p( ^> k-Y k 


l n =0 
\ 1 




n. ! I 
1 / 


R_1 M 1 I, M 1 U / N+N a 

V /"“l-vr-k n-l-k ' A n n , 

^ V t )P q E p, 1/n, ! 

k-0 k \nj-0 1 1 


N n n 
E 0 1 A /n,ll 

v° / 


(B. 17) 
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E(n^) 

P(o) 


R-l 

E 

k-0 


.M-l. k M-l-k 
<k )Pq 




. R M-l-R / 

)p q p, 


N-l 

E 

m«0 
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Therefore, the effective boundary is given, on average, by: 

M E(n.) 

E(Y) - -J- 1 - (B.19) 

Similarily, the second moment of the effective boundary (or the ntnnber of 
busy CS slots) is, after some calculation given by 


,, P(o) 
E(Y“) = — =“ 
F 


E ( M ^ 1 )p k q M_1 ' k pj[ A(S+K a -2)-K> 1 A(N+K a -1) 
k~0 


+(M R 1)pRqM_1 Pi A < N_2 > + P^N-l) 
m 

where A(m) ■ I P?Vi! 

i=0 


(B.20) 


(B.21) 


2 

The parameter E(Y ) is needed to evaluate the packet time delay using the 
model of a single server with "breakdowns" (see Eq.(2.14)). 
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4. OPTIMAL CIRCUIT SCHEDULING IN SSMB SYSTEMS: PERFORMANCE BOUNDS AND 
COMPUTATIONAL COMPLEXITY 

Arbitrary traffic matrices 
Interference (overlapping beam) constraints 
Multi-channel beams 
Different uplink/down configurations 
As indicated in earlier chapters, when system structure and constraints 
become more complex, the problems of design and performance analysis 
increase commensurately. Thus, we have generally limited our previous 
discussion of system designs to suboptimal ones when dealing with packet 
switched systems having such complexities as skewed traffic matrices and 
overlapping beams. Furthermore, our performance evaluation in these cases 
has been essentially by simulation rather than analysis, since in a packet- 
switching or integrated networks environment, any analytical models would be 
completely intractable. However, if we limit our attention to systems in 
which capacity is assigned on a slowly varying basis, as would be the case 
for dedicated or switched circuit assignments, it is possible to expand 
considerably the range of systems which are amenable to optimal design and 
analysis. In this chapter we shall examine the SS/TDMA scheduling problem in 
this context. The thrust of our work on this problem has been: 

1) To determine the computational complexity of the scheduling problems 
with a given set of constraints. 

2) To develop optional scheduling algorithms in cases where that is 
feasible, or suboptimal heuristics in cases where optimization is 


impractical 
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We consider system with arbitrary (i.s. skewed) traffic matrioes, with 
possible overlap (i.e. , interference) between beam and with one or aore of 
the following additional types of complexity: 

• different uplink/downlink beaa configurations 

- different power allocations in different beam 

- frequency reuse within beam by polarization 

- channelization within beam by PIMA, CDMA or other orthogonal 

multiplexing aohems 

- onboard demod/renod 

• fewer transponders than beams 

In the following sections, we givti a nummary of our results and those of 
others on these problems. The discussion emphasizes general computational 
features of the scheduling problem. The specifies of the algorithm can be 
found in the cited references. 

4.1. OPTIMAL SWITCHING FOR SYSTEMS WITH SINGLE CHANNEL BEAMS 

We consider a system operatng with a fixed TDMA frame structure in 
which slots are to be assigned on a permanent (or at least relatively long 
term) basis to traffic defined by a given traffic demand matrix D. The 
assignment of a single slot is assumed to represent the allocation of "one 
unit" of capacity. The "units" might be in the form of digitized voice 
channels (i.e. 64 kbps ohannels) or in any other convenient size. The entry 
djj (an integer) in the D matrix indicates the number of units of capacity 
required for transmissions from zone i to zone j. This capacity might be 
allocated on a permanent basis (if we are dealing with leased channels) or it 
might represent a pool of resources set aside for demand assignment in a 
circuit-switched manner. In either case the slot allocation problem is the 
same. We consider first, a system which has N uplinks, switched (through an 
r-f switch in the satellite) to N downlinks, so that the demand matrix <s Nx 
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N. (We treat below the oase where there may be a different number of uplinks 
and downlinks.) The optimal switching problem oonslsts of aooomaodatlng the 
specified traffic demand within a TDMA frame of ainiaum duration subject to 
the system constraints. These constraints inolude the possible conflicts 
within the switch, limitations on the number of active transponders, 
interference due to overlapping beams, and possibly others. A minimum 
duration frame corresponds to maximum transponder utilization. Thus, an 
optimized assignment corresponds to maximizing the throughput per frame. 

(In practice, the time duration of the frame might be fixed, in which case 
solving the optimization problem as stated here corresponds to packing a 
given demand as "tightly" as possible into the frame, thereby allowing as 
much as possible for additional demand, or more generally, freeing as much 
transponder capacity as possible for other use.) 

Assuming that there are N transponders available and that there are no 
other constraints except those of the switch, it can be shown that a frame of 
length 


where 


L > Max [max R if max Cj ] 


R i s j d ij 


V \ d ij 


<«.» 


is necessary to accommodate all traffic. For example, for the 4x4 
matrix of Fig. 4.1 the minimum possible frame size is L s 4. It can also be 
shown that this minimum length is achievable, and an algorithm has been 
developed for determining an optimal schedule. As an illustration, the 


matrix of Pig. 4.1 can bo scheduled as shown in Pig. 4.2. In tha figure, aaeh 
row aight ba oonsidered to represent tha transmission* of a glvan tranapondar 
assooiatad with that uplink. Note that oartain tranapondara ara inaotiva 
during part of tha frame (shown shaded in the figure), indicating that 100# 
transponder utilization is not (and cannot be) attained in this example. (Of 
course, the assignment is nevertheless optimal.) For the details of the 
optimal scheduling algorithm the reader ia referred to [iNU 79] . 

It is interesting to note that the assignment problem just described is 
related to a well-known problem in graph theory, known as the edge coloring 
problem in bipartite graphs. As will be seen below, many of the other 
problems we consider, involving more realistic constraints also have their 
counterparts in graph theory, suggesting that there is a l<.**ge body of 
literature to draw on in developing algorithms for scheduling and capacity 
allocation. 

The above problem can be generalized to the case where there are less 
active transponders than zones and the demand matrix is not square. A non- 
square matrix corresponds to a system with a different number of uplink and 
downlink zones. Such a configuration might result from a system design 
adapted to highly non-symmetrio traffic distributions, in whioh antenna gains 
and power allocations are adjusted to define geographic coverages differently 
for uplink and downlink transmissions. Or it might be a design especially 
adapted to certain special broadcast (multi-destination) requirements. A 
typical example is shown in Fig. 4.3. In this context, the optimal 
scheduling problem is stated as follows. Given a (generally rectangular) 
demand matrix D, find a schedule which assigns all traffio in a frame of 
minimum duration using K transponders. (The limitation to K transponders 
corresponds to a constraint that no more than K simultaneous transmission can 
occur during any slot.) 
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It oan be shown that tha minimum frame length neoessary to sohedtle all 
trafflo In tills oasa is given by 

L * Max [max R t , max C^, f|]] («l.2) 

where T a I d 

l.j lJ 

and f*x"] a least integer ,> x 

An algorithm for performing the scheduling with this value of L has been 
developed. As an example, with K a 2 transponders, the demand matrix of 
Fig. 4.4 can be scheduled within a frame of length L a 12, as shown in the 
figure. (This is the minimum possible frame length.) The details of the 
optimal scheduling algorithm in this oase may be found in [B0N801 . 

Sometimes it is advantageous to attempt to reduce to a minimum the 
number of times the switch is reconfigured within a frame. This might be the 
case for example, when switching requires a significant amount of time or 
power. Previous work [IT0771 has treated the scheduling problem with a more 
or less "informal” attempt to reduce the number of swltohings Ns. However, 
we have found that a true optimal soheduling approach in whioh a constraint 
on Ns is included in the formulation, generally turns out to be an 
exceedingly difficult combinatorial problem. We comment briefly here on the 
computational complexity of this problem as determined from our recent 
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studies, restricting ourselves to the ease of an M x M aatrlx D with K 
transponders. It is easily Mown that the nlnlnun number of switohings *min 
required for any D is 4 N. Thus one approaoh would be to determine **min 
and then attempt to determine an optimal sohedule subject to the constraint, 
N # * N oin* Son ® P* rtial results have appeared in the literature on this 
problem. [CAL801 We have determined that the problem is NP-c' .opiate, i.e. , 
its computational oomplexity in all likelihood grows exponentially with N. 

On the other hand, the teohniques mentioned above, for optimal soneduling 
without constraints on N are guaranteed to produoe N.*N 2 . These teohniques 
have a complexity of the order of N^. For constraints on Ns between these 
two extremes, (N and N ) it is uncertain at this point what the complexity of 
the problem is. This problem is still under study, and various heuristios 
have been developed for obtaining "good” schedules subject to switching 
constraints. 

4.2 MULTI-CHANNEL AND DUAL POLARIZED BEAMS 

So far it has been assumed that all uplink and downlink beams are 
of equal capacity, each permitting one unit of traffic to be transmitted in 
each time slot. In reality, eaoh beam is likely to be broken down into 
subchannels in various ways. Por example, Independent transmissions eould be 
superimposed in the same frequency band using two orthogonal polarizations. 

Or any one of a number of standard orthogonal multiplexing teohniques such as 
FDMA, CDMA or other analog and digital methods could be used to to create 
several independent channels within the same beam, the number of channels 
being a function of bandwidth, antenna gain and power allocated to the beam. 
Such channelization might be useful for oases in which heterogeneous 



I 

populations of Mali and largo stations exist within a s ingle boss, and/or ) 

the trafflo natrlx la highly skewed and asyMotrlo. In the lattar oase, 

beans with a large nuaber of ohannels eould be oonflgured to oover sones of 

high trafflo denand, while zones of lower demand would be served by beans 

with a lower nuaber of ohannels, A general nodal applloable to a wide range 

of suoh systems Is shown in Pig, 4,5, As shown in the figure, the i-th 

uplink is assuaed to be channelized Into subohannels eaoh of one unit 

capacity, and the J-th downlink channelized into subohannels of me unit 

» 

eaoh. The data streaas arriving at the satellite are deaultlplexed, so that i 

f 

I 

a total of subchannels are available at the input of the first switch, ’ 

to be connected to the set of Zcj subchannels on the downlinks. The 
configuration shown in the figure (Just one of many possible ones) uses a 

I 

single transponder for each subchannel. Thus the maximum number of ) \ 

subchannels that can be active in any one time slot will be K, the number of 

available transponders. More generally, K represents the maximum number of 

simultaneous transmissions of which the satellite is capable. This might be 

a function of switch, power or other limitations, in addition to limitations 

on the nuaber of transponders. He consider the general case of rectangular 

systems (M uplinks and N downlinks). Note that the demultiplexing and 

multiplexing night consist only of filtering and frequency translation (as in 

an FDMA system), or it might involve onboard demodulation and remodulation in 

more complex system. 

The problem of optimal scheduling in this systen is similar to that 
treated above Involving systems with a single channel per beam. (Illustrated 
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in Pig. 4.4). Tils time, however, the switch constraints are more coop lex. 

Now In a given time slot, a aaxiaua of g^ units of traffio oan be 
transait ced on uplink i, a aaxiaua of <*j units on downlink J and the 
total .ranamissions in any slot oannot exceed K. As before, we start with a 
(reotangular) demand matrix D with non-negative integer entries, and wish to 
find a schedule for this traffic in a minimum frame duration L, respecting 
the switch constraints mentioned above. It has been shown that a lower bound 


on frame duration for this 

oase is give, by 


i 

i 

| 

L * Max [" [ ~ 1 

. maxf\~] , 

oax ['—111 

\ 

(4,3) \ 

M k 1 

i 1*1 * 

J ‘vJ 


Furthermore an algorithm has been developed for constructing a schedule 
which achieves this bound. The reader is referred to [GOP 811 for details 
of the algorithm. He will oonfine ourselves here to a simple illustrative 
example. Fig. 4.6 shows a demand matrix for a three beam system having 
various channelizations of its up and down links as indicated by the 
S' s and c's in the figure, and having K « 4 transponders. For this 
example, the lower bound of (4.3) yields L > 5. By applying the algorithm a 
schedule can be constructed which achieves this bound, as shown in the 
figure. The traffic matrix D is broken down into a "feasible decomposition", 
which is a sum of matrices D^, eaoh of which represents the traffic 
transmitted in one tim* slot. The matrix in figure shows that in slot 1, 
a total of four units of traffio are transmitted (consistent with the 
constraint K « 4). It is distributed as follows: One subohannel in beam 1 
is transmitting back down to beam 1, and the other is connected to a 
subohannel in downlink 2. Simultaneously, the two subchannels in uplink 2 
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arc oonneoted back down to two subchannels in beam 2. This oompletely 
exhausts all subohannals in tha and downlines on beams 1 and 2, and all 
transponders, so that uplink 3 does not transnit any traffio in this slot. 

Tha remainder of tha schedule is defined in a similar fashion by matrices 02 * 
D 3 , Oj, and D^. 

4.3 ADJACENT BEAM INTERFERENCE 

In this section we confine attention to N x N systems. In 
general, if a large contiguous area is to be covered by a number of beams in 
a multibeam system, there will always be some Interference in the regions of 
zones affected by adjacent beam overlap, as illustrated by the shaded areas 
in Fig. 4.7. More generally, depending upon how well the antenna side lobes 
are suppressed, this interference may even spill over into non-adjacent 
zones. It is convenient to represent these interference effects by an 
"interference graph", the edges of which signify mutual interference between 
the vertices (zones) at their extremities. In a graph-theoretic sense such 
vertices are "adjacent", even though the corresponding zones may not be 
geographically adjacent. In the following, we will define "adjacent zone” 
interference in this more general sense. Fig. 4.8 shows an interference 
graph superimposed on the geographical interference picture shown in Fig. 

4.7. Note that in the case where interference only occurs between 
geographically adjacent zones, the interference graph is planar. As will be 
pointed out below, planar interference patterns are more easily handled than 
non-planar ones. 

It is dearly not feasible to simultaneously transmit to, or receive 
from two zones which mutually interfere. There are a number of ways of 
"rear- * “*ng" such interference. If transmissions to/' from all zones are 
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ass used to be reusing a common uplink and common downlink 
channel (as would be the ease, for example. In an Ixl system reusing the 
sane frequency band N times, without exploiting orthogonal polarizations or 
other methods of channelization within a bean) then interference can only be 
resolved by scheduling the TDMA frame so that at most one of a mutually 
Interfering set of beams is transmitting (receiving) at a time. In the 
example of Fig. 4.8, zones 2 and 3 can be scheduled simultaneously, but if 1 
or 4 Is active, all others must be inactive. Clearly, a schedule for a given 
demand matrix in a system which requires resolution of adjacent zone 
interference will yield a longer frame length than a schedule for a non- 
interfering system. 

A second way of resolving interference is to make use of several 
orthogonal channels of the type proposed in Sec. 4.2. For example, if two 
orthogonal polarizations are available, pairs of interfering zones could be 
simultaneously active provided that they were using different polarizations. 
Similarly, if two or more separate frequency bands (or other forms of 
orthogonal channels) were available, they could be used for simultaneous 
transmissions to or from a mutually interfering set of beams. In general, 
then, we resolve interference by requiring that transmissions to or from a 
set of interfering beams be orthogonal. If orthogonality in time is used, as 
described above, then the time frame generally must be expanded to resolve 
interference. If orthogonality in frequency is used we are thereby "cutting 
up" the available frequency spectrum and limiting its reuse. In either case, 
resolution of interference will necessarily reduce the efficiency of 
utilization of the available communication resources. Since efficient 
resource allocation is a prime consideration in SS/TDMA systems, we have 
studied possible algorithms for optimal scheduling in the face of adjacent 


zone interference. 
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Consider, first, the problem of resolving int erferen c e by the use of two 
polerizetions. The overall problee is: given e traffic eatruc D, schedule 
this traffic in einieue tine subject to both the switch/transponder 
constraints and the interference constraints. 

One way to approach this problee is via a two-step approach: 

1) Taking into account the traffic eatrix D, assign polarizations to the 
various transmissions which reduce the interference conflicts to a einieue. 

2) If all interference has not be resolved in step (1), expand the 
resulting schedule to a larger ties frame by adding time slots in order to 
separate in time, the remaining interfering traffic. 

In the case where the interference graph is planar we have formulated an 
algorithm for step (1) which requires polynomial time; i.e. , it requires a 
computational effort which grows as a power of N. The algorithm is based on 
a procedure of Hadlock for solving a problem in graph theory known as the 
Max-cut problem. 

In the case of non -planar interference graphs, and/or in cases where 
more than two orthogonal channels are available for resolving interference, 
we have shown that the scheduling problem is NP -complete. Since it would be 
impractical to implement an optimal scheduling algorithm in this more general 
case, we have developed several heuristics for constructing subop tlmal 
schedules. By establishing a lower bound on the frame length for the optimal 
schedule with interference constraints, we have been able to estimate the 
efficiencies of these heuristics. Our simulation results averaged over a 
fairly large number of tests show that the frame lengths based mi heuristic 
approaches average only about 15 % more than the (not necessarily achievable) 
lower bound. This work is still in progress. 
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Fig. 4.8 Interference graph. 
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